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1. Introduction 

These lectures present an introduction to quantum cosmology for an audience 
consisting for a large part of astronomers, and also a number of particle physicists. 
As such, the material covered will invariably overlap with that of similar intro- 
ductory reviews 1 ' 2 , although I hope to emphasise, where possible, those aspects of 
quantum cosmology which are of most interest to astronomers. 

This is not an easy task - quantum cosmology is not often discussed at Summer 
Schools such as this where there is a large emphasis on astrophysical phenomenology, 
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for the very good reason that the ideas involved are at present rather tentative, and 
quantitative predictions are thus much more difficult to arrive at. Nevertheless 
we should remember that not too long ago the whole enterprise of cosmology was 
viewed as the realm of wild speculation. Vladimir Lukash remarked that when he 
started out in research the advice he was given was that "cosmology was all right 
for someone like Zeldovich to potter around with after he had already established 
himself by producing a body of serious work, but it was inappropriate for a young 
physicist embarking on a career". Happily the situation is quite different today - 
I am sure the earlier lectures in this School will have convinced you that modern 
cosmology is a hardcore quantitative science, and that with the new technology 
and techniques now being developed we can expect to accurately measure all the 
important cosmological parameters within the next decade and thus enter into a 
"Golden Era" for cosmological research. 

Quantum cosmology, however, still enjoys the sort of status that all of cosmology 
had until not so very long ago: essentially it is a dangerous field to work in if you 
hope to get a permanent job. I hope to convince you nevertheless that quantum 
cosmology represents a vitally important frontier of research, and that although it is 
by nature somewhat speculative, such speculations are vital if we are to understand 
the entire history of the universe. 

On the face of it the very words "quantum" and "cosmology" do appear to some 
physicists to be inherently incompatible. We usually think of cosmology in terms of 
the very large scale structure of the universe, and of quantum phenomena in terms 
of the very small. However, if the hot big bang is the correct description of the 
universe - which we can safely assume given the overwhelming evidence described 
in the earlier lectures - then the universe did start out incredibly small, and there 
must have been an epoch when quantum mechanics applied to the universe as a 
whole. 

There are people who would take issue with this. In the standard Copenhagen 
interpretation of quantum mechanics one always has a classical world in which the 
quantum one is embedded. We have observers who make measurements - the ob- 
servers themselves are well described by classical physics. If the whole universe is to 
be treated as a quantum system one does not have such a luxury, and some would 
argue that our conventional ideas about quantum physics cease to make sense. Yet 
if quantum mechanics is a universal theory then it is inevitable that some form of 
"quantum cosmology" was important at the earliest of conceivable times, namely 
the Planck time, t Planck = (Gh) 1/2 /c 5 / 2 = 5.4 x 10" 44 sec, (equivalent to 10 19 GeV 
as an energy, or 1.6 x 10~ 35 m as a length). At such scales, where the Compton 
wavelength of a particle is roughly equal to its gravitational (Schwarzschild) radius, 
irreducible quantum fluctuations render the classical concept of spacetime mean- 
ingless. Whether or not our current efforts at constructing a theory of quantum 
cosmology are physically valid is therefore really a question of whether our current 
understanding of quantum physics is adequate for considering the description of 
processes at the very beginning of the universe, or whether quantum mechanics 



2 



itself has to be revised at some level. Such a question can really only be answered 
by extensive work on the problem. 

Setting aside the question of the fundamentals of quantum mechanics, let us 
briefly review the problems which are left unanswered in the standard hot big bang 
scenario. These are: 

1. The value ofn = N, IN , . : the exact value of this parameter is unexplained 

J ' baryons' photons A A 

in the hot big bang but is crucial in determining abundances of light elements 
through primordial nucleosynthesis. 

2. The horizon problem: the isotropy of the cosmic microwave background radiation 
indicates that all regions of the sky must have been in thermal contact at some time 
in the past. However, in the standard Friedmann-Robertson- Walker (FRW) models 
regions separated by more than a couple of degrees have non-intersecting particle 
horizons - i.e., they cannot have been in causal (and hence thermal) contact. 

3. The flatness problem: given the range of possible values of the ratio of the density 
of the universe to the critical density at the present epoch, Q , then the FRW models 
predict that Q(t) must have been incredibly close to the spatially flat case Q(t) ~ 1 
at early times; e.g., assuming fi Q > 0.1 we find \ft — 1| < 10~ 26 at the lepton era, 
and \ ft - 1| < 10" 53 at the GUT era. 

4. The unwanted relic problem: models of the early universe which involve phase 
transitions often produce copious amounts of topological defects, such as monopoles 
produced at the GUT scale. If one puts the numbers in one finds that the density of 
such relics is so great that they would exceed the critical density by such a margin 
in the standard FRW models that the universe should have ended long ago! 

5. The origin of density perturbations is unexplained. 

6. The arrow of time is a physical mystery. On the one hand the laws of physics 
are CPT-invariant, and on the other there is a thermodynamic arrow of time, as 
prescribed by the Second Law of Thermodynamics, and it appears to match the 
cosmological arrow of time, as prescribed by the expansion of the universe. 

7. The initial conditions of the universe must be put in by hand, rather than being 
physically prescribed. 

The first four problems on the list are ones that can be explained without 
appealing to quantum cosmology. The value of n (problem 1) is predicted by models 
of baryogenesis, which typically take place at the GUT scale. Problems 2-4 are 
solved by the inflationary universe scenario: an early phase of exponential expansion 
of the universe drastically changes the past light cone, thereby removing the horizon 
problem, while also driving Q close to unity, and diluting unwanted relics to such 
very low densities that they are close to unobservable. 

Problems 5-7, on the other hand, are of a nature which is beyond the scope 
of the inflationary universe scenario to satisfactorily explain. Inflation provides a 
mechanism whereby initial small "quantum" perturbations are inflated to all length 
scales to form a scale-free Harrison-Zeldovich spectrum, but it does not address 
the question as to exactly how these perturbations arise. Furthermore, a typical 
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model of the very early universe might possess both inflationary and non-inflationary 
solutions, so that the precise initial conditions of the universe can be crucial for 
determining whether the universe undergoes a period of inflation sufficiently long to 
be consistent with observation. The length of the period of inflation is precisely the 
sort of quantitative result that we might hope quantum cosmology should provide. 
Questions such as the origin of the arrow of time might appear to be of a more 
philosophical nature - however, quantum cosmology should provide a calculational 
framework in which such questions can begin to be addressed. 

1.1. Quantum cosmology and quantum gravity 

Quantum cosmology is perhaps most properly viewed as one attempt among 
many to grapple with the question of finding a quantum theory of gravity. As a 
field theory general relativity is not perturbatively renormalisable, and attempts to 
reconcile general relativity with quantum physics have not yet succeeded despite 
the attentions of at least one generation of physicists. It is perhaps not surprising 
that the problem is such a difficult one since general relativity is a theory about the 
large scale structure of spacetime, and to quantise it we have to quantise spacetime 
itself rather than simply quantising fields that live in spacetime. 

Many ideas have been considered in the quest for a fundamental quantum the- 
ory of gravity - whether or not these ideas have brought us closer to that goal is 
difficult to say without the benefit of hindsight. However, such ideas have certainly 
profoundly increased our knowledge about the nature of possible physical theories. 
Some important areas of research have included: 

1. Supergravity 3 ^. Using supersymmetry, a symmetry between fermions and 
bosons, one can enlarge the gravitational degrees of freedom to include one or more 
spin-| gravitinos, in addition to the spin-2 graviton, g^ v . Such a symmetry can 
cure some but not all of the divergences of perturbative quantum general relativity*. 
In particular, while pure Einstein gravity is perturbatively non-renormalisable at 
two loops 7 , or at one loop if interacting with matter 8-10 , in the case of supergravity 
renormalisability fails only at the 3- loop level 11 ^. 

2. Superstring theory 13 . Progress can be made if in addition to using supersymme- 
try, one constructs a theory in which the fundamental objects have an extension 
rather than being point-like: a theory of strings rather than particles. Much inter- 
est in string theory was generated in the mid 1980s with the discovery that certain 
string theories appear to be finite at each order of perturbation theory. In some 

* For details of the application of perturbative techniques to quantum cosmology, both with and without 
supersymmetry, see [6]. 

^ The status of the result concerning the 3-loop divergence of supergravity is not quite as rigorous as the 
other examples mentioned, as the complete 3-loop calculation has not been done. However, it is known that 
a 3-loop counterterm exists for all extended supergravities and there is no reason to expect the coefficient 
of the counterterm to be zero, making 3-loop finiteness extremely unlikely. For a review of the ultraviolet 
properties of supersymmetric field theories see [12]. 
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sense stringy physics "smears out" the problems associated with pointlike interac- 
tions. The entire sum of all terms in the perturbation expansion diverges in the 
case of the bosonic string 14 , however, and it is believed that similar results should 
apply to the superstring®. Furthermore, despite what many see as the mathemati- 
cal beauty of string theory, there has unfortunately as yet been no definitive success 
in deriving concrete phenomenological predictions. 

3. Non-perturbative canonical quantum gravity 16,17 . The fact that general relativity 
is not perturbatively renormalisable might simply be a failure of fiat spacetime 
quantum field theory techniques to deal with such an inherently non-linear theory, 
rather than reflecting an inherent incompatibility between general relativity and 
quantum physics. Given the divergence of the string perturbation series mentioned 
above, a non-perturbative formulation of string theory would also be desirable. As a 
starting point a systematic investigation of a non-perturbative canonical formalism 
based on general relativity could provide deep insights into quantum gravity. Such a 
programme has been investigated by Ashtekar and coworkers, mainly from the mid 
1980s onwards. One principal difference from the canonical formalism which I will 
describe in section 2 is that instead of taking the metric to be the fundamental object 
to describe the quantum dynamics, one bases such a dynamics on a connection, in 
this case an SL(2, C) spin connection 18 ' 19 . If one integrates this connection around 
a closed loop one arrives at "loop variables", which might be considered to be 
analogous to the magnetic flux, ® 7 obtained by integrating the electromagnetic 
gauge potential, A^, around a closed loop. In the "loop representation" quantum 
states are represented by functionals of such loops on a 3-manifold 20 ' 21 , rather 
than by functionals of classical fields. Although a number of technical difficulties 
remain, considerable progress has been made with the Ashtekar formulation, and 
the reformulation of quantum cosmology in the Ashtekar framework 22-25 could be 
an area for some interesting future work. 

4. Alternative models of spacetime. All the above approaches assume that the basic 
quantum variables, be they a metric or a connection, are defined on differentiable 
manifolds. Given that it is highly possible that "something strange" happens at the 
Planck length it is plausible that one might have to abandon this assumption in order 
to effectively describe the quantum dynamics of gravity. A number of ideas have 
been considered on these lines. One framework which has been widely used both 
for numerical relativity and studies of quantum gravity is that of Regge calculus 26 lA 
in which one replaces smooth manifolds by spaces consisting of piecewise linear 
simplicial blocks. Naturally, other possibilities for discretised spacetime structure 
also exist - causal sets 28 being another example which has not yet been so widely 
explored. A further possibility is topological quantisation whereby one replaces a 
manifold by a set and quantises all topologies on that set 29 ' 30 . 

(5) 

The question of finiteness of superstring perturbation theory is a difficult technical question, which 

has still to be resolved - see, e.g., [15]. 
A For a brief review and extensive bibliography see [27]. 
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Many of the above alternatives fall into the category of being attempts to con- 
struct a fundamental theory of quantum gravity. The current quantum cosmology 
programme is not quite as ambitious. One begins by making the assumption that 
whatever the exact nature of the fundamental theory of quantum gravity is, in its 
semiclassical limit it should agree with the semiclassical limit of a canonical quan- 
tum formalism based on general relativity alone. Thus we study the semiclassical 
properties of quantum gravity based solely on Einstein's theory, or some suitable 
modification of it. 

Clearly, any predictions made from such a foundation must be treated with 
caution. In particular, a new fundamental theory of quantum gravity might in- 
troduce radically new physical processes at an energy scale relevant for cosmology. 
String theory, for example, introduces its own fundamental scale which expressed 
as a temperature (the Hagedorn temperature) is given by: 



the constant a' being the Regge slope parameter, which is inversely proportional 
to the string tension. It is not known exactly what the value of ot is; however, 
if 7" Hagodorn is comparable to or significantly lower than the Planck scale, then it 
is clear that fundamentally "stringy" processes will be very important in the very 
early universe, if string theory is indeed the "ultimate" theory. 

Although new fundamental physics could drastically change the predictions of 
quantum cosmology, I believe nevertheless that studying quantum cosmology based 
even just on Einstein's theory is an important activity. General relativity is a 
remarkably successful theory; in seeking to replace it by something better it is im- 
portant that we study processes at the limit of its applicability, thereby challenging 
our understanding. General relativity is limited by the Planck scale - the physical 
arenas in which this scale is approached include: (i) very small black holes; (ii) the 
very early universe. Since it seems unlikely that we will ever be able to create ener- 
gies of order 10 19 GeV in the laboratory, the consequences of quantum gravity for 
the physics of the very early universe will remain the one way of indirectly "testing" 
it, at least for the foreseeable future. 

It is thus important that we consider quantum gravity in a cosmo logical context. 
Even if our current attempts do not fully reach the mark, in that we do not yet 
have a fully-fledged quantum theory of gravity, they nonetheless constitute a vital 
part of the process of trying to find such a theory. 

1.2. A brief history of quantum cosmology 

The quantum cosmology programme which I will describe in these lectures has 
gone through three main identifiable phases to date: 



Hagedorn 




(1.1) 
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1. Defining the problem. The canonical formalism, including the definition of the 
wavefunction of the universe, its configuration space - superspace - and its evo- 
lution according to the Wheeler-DeWitt equation, was set up in the late 1960s 31-35 . 

2. Boundary conditions. Quantum cosmology research went into something of a lull 
during the 1970s but was revived in the mid 1980s when the question of putting 
appropriate boundary conditions on the wavefunction of the universe was treated 
seriously. The idea is that such boundary conditions should describe the "creation 
of the universe from nothing" 36 ' 37 , where nothing means the absence of space and 
time. A number of proposals for such boundary conditions emerged - two major 
contenders being the "no-boundary" proposal of Hartle and Hawking 38 ' 39 and the 
"tunneling" proposal advocated by Vilenkin 40-42 . 

3. Quantum decoherence. The mechanism of the transition from quantum physics to 
classical physics ("quantum decoherence") becomes vitally important when quan- 
tum physics is applied to the universe as a whole. The issues involved have begun 
to occupy many researchers in the early 1990s T . 

Two other important areas of quantum cosmology (or related) research have been: 
(i) quantum wormholes and "baby universes"; (ii) supersymmetric quantum cos- 
mology. 

Quantum wormholes* were extremely fashionable in the particle physics com- 
munity in the years 1988-1990. Such states arise when one considers topology 
change in the path integral formulation of quantum gravity: quantum wormholes 
are instanton solutions which play an important role in the Euclidean path inte- 
gral. One deals directly with a "third quantised" formalism, (i.e., quantum field 
theory over superspace), which includes operators that create and destroy universes 
- so-called baby universes. Much of the excitement in the late 1980s was associated 
with the idea that such processes could fix the fundamental constants of nature - 
in particular, driving the cosmological constant to zero. 

Supersymmetric quantum cosmology has emerged recently as one of the most ac- 
tive areas of current research^. In considering the quantum creation of the universe 
we are of course dealing with the very earliest epochs of the universe's existence, at 
which time it is believed that supersymmetry would not yet be broken. The inclu- 
sion of supersymmetry could therefore be vital from the point of view of physical 
consistency. 

Since the focus of this School is on cosmology, my intention in these lectures 
is to cover topics 1 and 2 above, and then to proceed to discuss the predictions of 
quantum cosmology. The third topic of quantum decoherence raises questions which 
have not been resolved even in ordinary quantum mechanics, since the question of 
decoherence really amounts to understanding what happens during the "collapse 
of the wavefunction". Although this is a fascinating issue it has more to do with 

T See [43] for a review and [44] for a collection of recent papers on the subject. 

* See [45] for a review. 

^ See [46] and references therein. 
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the fundamentals of quantum mechanics than directly with cosmology. Likewise, 
I will only briefly touch upon quantum wormholes and supersymmetric quantum 
cosmology, as these still in their infancy, and one is still at the point of 

trying to resolve basic questions concerning quantum gravity. I hope the reader will 
not be disappointed by this - however, given the vast scope of quantum gravity and 
quantum cosmology one must necessarily be rather selective. 



2. Hamiltonian Formulation of General Relativity 

2.1. The 3 + 1 decomposition 

In order to discuss quantum cosmology a fair amount of technical machinery 
is required. In the canonical formulation we begin by making a 3 + 1-split of the 
4-dimensional spacetime manifold, M. , which will describe the universe, foliating it 
into spatial hypersurfaces, Ef, labeled by a global time function, t. Thus we take 
the 4-dimensional metric to be given by* 



ds 2 = g liV dx' i dx v = -uj° <g> + hijU* <8> 



(2.1) 



where 



uj° = Afdt 

<J = dx'+Af^t. 



(2.2) 



Such a decomposition is possible in general if the manifold M is globally hyperbolic^ 7 . 
The quantity J\f(t, x k ) is called the lapse function - it measures the difference be- 
tween the coordinate time, t, and proper time, r, on curves normal to the hyper- 
surfaces E t , the normal being n a = (—A/", 0,0,0) in the above coordinates. The 
quantity M l (t, x k ) is called the shift vector - it measures the difference between a 
spatial point, p, and the point one would reach if instead of following p from one 
hypersurface to the next one followed a curve tangent to the normal n. That is 
to say, the spatial coordinates are "comoving" if Af l = 0. Finally, hij{t,x k ) is the 
intrinsic metric (or first fundamental form) induced on the spatial hypersurfaces by 
the full 4-dimensional metric, g^ v . In components we have 



(9^) = 



-N 2 +N k N k jsfj 
Mi h 



•J 



(2.3) 



with inverse 



or ) = 



JT 1 



Af 2 

t,ij Af i Af j 



(2.4) 



where h tJ is the inverse to hij, and the intrinsic metric is used to lower and raise 



spatial indices: A/^A/fc = h jk UjUk = h jk N' j N' k etc. 



Wc use a ( — h ++) Lorcntzian metric signature, and natural units in which c = h = 1. 
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x l +dx l 

Fig. 1: The 3 + 1 decomposition of the manifold, with lapse function, N, and shift 
vector, N l . 

One can construct an intrinsic curvature tensor 3 R l j kl (h) from the intrinsic 
metric alone - this of course describes the curvature intrinsic to the hypersurfaces 
Sf. One can also define an extrinsic curvature, (or second fundamental form) , which 
describes how the spatial hypersurfaces curve with respect to the 4-dimensional 
spacetime manifold within which they are embedded. This is given by 

Kij = — ni-j = — r ^n 

where a semicolon denotes covariant differentiation with respect to the 4- metric, g^, 
and a vertical bar denotes covariant differentiation with respect to the 3-metric, hij: 
= Mi,, -T k i5 M k etc. 

For a given foliation of M. by spatial hypersurfaces, S t , it is always possible to 
choose Gaussian normal coordinates, in which 

ds 2 = -dt 2 + h XJ dx l dx j . (2.6) 

These are comoving coordinates (Af l = 0) with the additional property that t is 
the proper time parameter (A/" = 1). This is the standard "gauge choice" that is 
made in classical cosmology, and in such coordinates = —hij, where dot denotes 



(2.5) 



differentiation with respect to t. In making the 3 + 1 decomposition, however, we 
are only free to make a specific choice of coordinates such as (2.6) after variation of 
the action if we want to be sure to obtain Einstein's equations, and thus we must 
retain the lapse and shift function for the time being. 



2.2. The action 



A relevant action for use in quantum cosmology is that of Einstein gravity plus 
a possible cosmological term, A, and matter, given by 



4k 2 



J d 4 x^ ( 4 R - 2A) + 2 J d 3 xVh~K 



M 



dM 



+ S, 



matter' 



(2.7) 



where k 2 = 4tvG = 47rmpi anck , K = K l i is the trace of the extrinsic curvature, 
and for many simple models the matter is specified by a single scalar field, with 
potential, V(<E>), 



S .. = I d Xy 
matter / 

M 



-yd^dv® - v($) 



(2- 



The boundary term 48 in (2.7) does not of course contribute to the classical field 
equations, and this term is usually omitted in a first course on general relativity. 
However, in quantum physics we are often interested in phenomena which occur 
when the classical field equations do not apply, that is "off-shell", and thus it is 
vitally important to retain the surface term. The simple matter action (2.8) given 
here should simply be seen as being representative of the type of matter action 
one might consider. Although the example given by (2.8) is sufficient for studying 
many inflationary universe models, many other alternatives might also be of interest, 
such as extra matter from a supergravity multiplet, or matter corresponding to the 
low-energy limit of string theory. In the latter case, if one works in the "string 
conformal frame" it is also necessary to alter the gravitational part of the action, 
as one characteristic of string theory is the presence of the scalar dilaton, which 
couples universally to matter (at least perturbatively) . In that case (2.7), (2.8) 
would be replaced by 



S = 



4k 2 



J d^v^e" 2 * ( 4 R + Ag^d^dv® - 8V($) + . . .) 

+ 2 j d 3 xVhe~ 2 *K 



dM 



(2.9) 



where the ellipsis denotes any additional matter degrees of freedom, and we have 
allowed for the possibility of the generation of a dilaton potential, V(<&), via some 
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non-perturbative symmetry breaking mechanism. However, for simplicity the spe- 
cific examples we will deal with here will be confined to models of the type (2.7), 
(2.8). 

We now wish to express (2.7), (2.8) in terms of the variables of the 3 + 1 split. 
One can show that 

*R = 3 R-2R nn + K 2 -K ij Kij, (2.10) 
where S R is the Ricci scalar of the intrinsic 3-geometry, and 

Run = R^n 13 = -K^Kn + K 2 + (n a K + a a ). a , (2.11) 

with a a = n^n^.p. Combining (2.7), (2.10) and (2.11), and noting that the bound- 
ary integral involving a a vanishes identically since n° l a cc = 0, we obtain 

S = JdtL = -^Jdt d'xAfVh (KyK** -K 2 + *R- 2A) + S mattc , (2.12) 

As in the Hamiltonian formulation of field theory we define canonical momenta in 
the standard fashion 

* s it = -5 - **) • <w 

^ = 0, (2,5) 

The fact that the momenta conjugate to M and Mi vanish means that we are dealing 
with primary constraints in Dirac's terminology 49 ' 50 . 

If we use A/", Mi, hij, <E> and their conjugate momenta as the basic variables we 
obtain^ a Hamiltonian 



H = J d 3 x (n°M + nWi + n ij hij + 7t*$) - L 
= J d 3 x (n°M + nWi + MH + NiH^j , 



(2.17) 



where 



^ Details of all missing steps in this section will be provided upon request in a plain brown envelope. 
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, (2.18) 



W = 



Vh 



2k 2 

= -27t y |j . + ^, i 7t., (2.19) 

and 

Qijki = \h~ 1 / 2 (hithji + huhjk — hijhki) (2.20) 

is the DeWitt metric 33 . The hats in (2.18) and (2.19) denote orthonormal frame 
components of the Einstein and energy-momentum tensors. In terms of these vari- 
ables the action (2.12) then becomes 

S = J dt d 3 x (n°Af + rtAfi -NH- AfiH^J . (2.21) 

If we vary (2.21) with respect to 7t lJ and 7t # we obtain their defining relations (2.13) 
and (2.14). The lapse and shift functions now act as Lagrange multipliers; variation 
of (2.21) with respect to the lapse function, A/", yields the Hamiltonian constraint 

H = 0, (2.22) 

while variation of (2.21) with respect to the shift vector, A/^, yields the momentum 
constraint 

W = 0. (2.23) 

From (2.18) and (2.19) it is clear that these constraints are simply the (00) and 
(Oi) parts of Einstein's equations. In Dirac's terminology these are secondary or 
dynamical constraints. 

The 3+1 decomposition of our spacetime looks to be very counterintuitive to the 
usual ideas of general relativity. This is so by choice. Einstein described the l.h.s. of 
his field equations which encodes the 4-dimensional geometry as a "hall of marble" , 
which does not encourage people to tamper with it. However, to quantise spacetime 
we must do just that: we must deconstruct spacetime and replace it by something 
else. Thus far we have not done that - we have simply split our 4-dimensional 
manifold into a sequence of spatial hypersurfaces, E t . Time is the natural variable 
to base this split upon since it plays a special role in quantum mechanics - it is a 
parameter rather than an operator. 

Classically the evolution of one spatial hypersurface to the next is completely 
well-defined, (provided that the manifold, A4, is globally hyperbolic), and given 
initial data, hij, <fr, on an initial hypersurface, E, we can use the Cauchy develop- 
ment to stitch the hypersurfaces E t together to recover the 4-dimensional manifold, 
M.. To quantise the theory, however, we want to perform a path integral over all 
geometries, not just the classically allowed ones. Thus we must consider sequences 
of geometries at the quantum level which cannot be stitched together in a regular 
Cauchy development to form a 4-manifold which solves Einstein's equations. (See 
Fig. 2.) We must therefore abandon Einstein's "hall of marble" - spacetime is no 
longer a fundamental object. 
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Fig. 2: Quantum geometrodynamics: in addition to the classical Cauchy devel- 
opment from S to £' (left), the path integral includes a sum over all 4- manifolds 
which interpolate between the initial and final configurations. The weighting by e lS 
means that the greater the number of classically forbidden 3-geometries (shaded 
slices) is in the interpolating 4-manifold, the smaller its contribution is to the path 
integral. 

As was mentioned in the Introduction the "deconstruction" of spacetime that 
we adopt here is probably the most conservative choice we could make. Even though 
we abandon the notion of spacetime in discussing the quantum dynamics of gravity, 
our fundamental objects are still defined on regular 3-manifolds, E. A more radical 
departure would be to replace these spatial hypersurfaces by some more general 
set by using, for example, the ideas of Regge calculus, causal sets or topological 
quantisation. 



3. Quantisation 

3.1. Superspace 

As a prelude to the canonical quantisation of gravity let us first introduce the 
relevant configuration space on which the quantum dynamics will be defined. 

Consider the space of all Riemannian 3-metrics and matter configurations on 
the spatial hypersurfaces, E, 

Riem (E) = {hy (x), | x G E} . (3.1) 

This is an infinite-dimensional space, on account of the label x = {x 1 }, which 
specifies the point on the hypersurface, but there are a finite number of degrees 
of freedom at each point, x G E. In fact, we are really interested in geometry 
here and configurations which can be related to each other by a diffeomorphism, 
i.e., a coordinate transformation, should be considered to be equivalent since their 
intrinsic geometry is the same. Thus we factor out by diffeomorphisms on the spatial 
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hypersurfaces and identify superspace* as 



Riem (E) 
Diffo(S) ' 

where the subscript zero denotes the fact that we consider only diffeomorphisms 
which are connected to the identity. This infinite-dimensional space will be the 
basic configuration space of quantum cosmology. 

The DeWitt metric (2.20) then provides a metric on superspace which we can 
write as 

G AB (x) = g {ijm) (x), (3.2) 

where the indices A, B run over the independent components of the intrinsic metric 

hij : 

A, B E {^11? ^12' ^135 ^22' ^23' ^33}' 

The DeWitt metric has a ( — h + + ++) signature at each point 16E, regardless 
of the signature of the spacetime metric, g^ u . To incorporate all the degrees of 
freedom, we also need to extend the range of the indices A, B to include the matter 
fields, by appropriately defining G^(x) (and other components if more than one 
matter field is present), thereby obtaining the full supermetric. 

An inverse DeWitt metric, Q AB (x) = G^ kl \x), can be defined by the require- 
ment 

G imk %kl)(mn) = \ ( V' n + W , (3-3) 

which gives 

g(ij)(kl) = i ^ ( h ik h ji + h u h jk _ 2h^h kl ) . (3.4) 



3.2. Canonical quantisation 

We will perform canonical quantisation by taking the wavefunction of the uni- 
verse, ty[hij,$>], to be a functional on superspace. Unlike the case of ordinary 
quantum mechanics, the wavefunction, does not depend explicitly on time here. 
This is related to the fact that general relativity is an "already parametrised" theory 
- the original Einstein-Hilbert action is time-reparametrisation invariant. Time is 
contained implicitly in the dynamical variables, and <E>. 

According to Dirac's quantisation procedure 49 we make the following replace- 
ments for the canonical momenta 

Shij 5$ 5Af 8Ni 



The use of the terminology "superspace" for the configuration space of quantum cosmology pre- 
dates the discovery of supersymmetry, and of course is completely different from the combined manifold of 
commuting and anticommuting coordinates which is called "superspace" in supersymmetric theories. 
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and then demand that the wavefunction is annihilated by the operator version of 
the constraints. For the primary constraints we have 



TtVjF = 



n l \& = -i 



5^ 



(3.6) 



0, 



which implies that \& is independent of Af and N" 1 . The momentum constraint yields 

5# 



H 1 ^ = 



Sh 



(3.7) 



Using (3.7) one can show 51 that \1/ is the same for configurations {hij (x), <&(#)} 
which are related by a coordinate transformation in the spatial hypersurface, E, 
which accords with the rationale for factoring out by diffeomorphisms in our defi- 
nition of superspace. Finally, the Hamiltonian constraint yields 



ijkl 



ShijShki ' 4k 2 V 



+ 



- a R + 2A + 4k z T 



2^00 



* = 0, 



where for our scalar field example 

-1 d 2 



H00 



(3.8) 



(3.9) 



Eq. (3.8) is known as the Wheeler-DeWitt equation 32,33 . In fact, it is not a single 
equation but is actually one equation at each point, x G E. It is a second-order 
hyperbolic functional differential equation on superspace. On account of factor- 
ordering ambiguities it is not completely well-defined, although there exist "natural" 
choices 33 ' 52 of ordering for which the derivative pieces become a Laplacian in the 
supermetric. 



3.3. Path integral quantisation 

An alternative to canonical quantisation which perhaps better accommodates 
an intuitive understanding (c.f. Fig. 2) of the quantisation procedure is the path 
integral approach. Path integral techniques in quantum gravity were pioneered in 
the late 1970s 53 ' 54 . The starting point for this is Feynman's idea that the amplitude 
to go from one state with an intrinsic metric, h^, and matter configuration, $, on an 
initial hypersurface, E, to another state with metric, h\ ■ and matter configuration 
$' on an final hypersurface, E', is given by a functional integral of e* s over all 4- 
geometries, g^, and matter configurations, cj), which interpolate between the initial 
and final configurations 

(ti ij ,&,v'\h ij ,$,v) = Y, [ DgDcK 5 ^-*]. (3.10) 
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In ordinary quantum field theory in flat spacetime the path integral suffers 
from the difficulty that since the action S[g^ v , 4>] is finite the integral oscillates and 
therefore fails to converge. Furthermore, the solution which extremises the action is 
given by solving a hyperbolic equation between initial and final boundary surfaces - 
which is not a mathematically well-posed problem, and may have either no solutions 
or an infinite number of solutions. To deal with this problem one performs a "Wick 
rotation" to "imaginary time" t — > —it and considers a path integral formulated 
in terms of the Euclidean action, / = — iS. The action is then positive-definite, 
so that the path integral is exponentially damped and should converge. Also the 
problem of finding the extremum becomes that of solving an elliptic equation with 
boundary conditions, and this is well-posed. 

One may attempt to apply a similar approach to quantum gravity, replacing S 
in (3.10) by the Euclidean action^ /[g^, <t>] = — iS[g^ v , c|>], and taking the sum in 
(3.10) to be over all metrics with signature (+ + +- 1-), which induce the appropriate 
3-metrics hij and h'^ on the past and future hypersurfaces. This approach to 
quantum gravity has had some important successes - most notably, it provides: 
(i) an elegant way of deriving the thermodynamic properties of black holes 55-57 ; 
and (ii) a natural means for discussing the effects of gravitational instantons 58-60 . 
Gravitational instantons have been found to provide the dominant contribution to 
the path integral in processes such as pair creation of charged black holes in a 
magnetic field 61-63 , and therefore provide a means of calculating the rates of such 
processes semiclassically. 

The problems associated with the Euclidean approach to quantum gravity are 
considerable, however. Firstly, unlike ordinary field theories* such as Yang-Mills 
theory the gravitational action is not positive-definite 64 , and thus the path integral 
does not converge if one restricts the sum to real Euclidean-signature metrics. To 
make the path integral converge it is necessary to include complex metrics in the 
sum ,65 . However, there is no unique contour to integrate along in superspace 66-68 
and the result one obtains for the path integral may depend crucially on the contour 
that is chosen 66 . Furthermore, the measure in (3.10) is ill-defined. It is really only 
in the last ten years that mathematicians have succeeded in making path integration 
in ordinary quantum field theory rigorously defined. Clearly, we may have to wait 
some time before the same can be said of path integrals in quantum gravity. 

The physicists' approach is to set aside the issues involved in making the for- 
malism rigorous and to see what can be learned nevertheless. We thus take the 

ES 

Strictly speaking one should call this the Riemannian action, since "Euclidean" spaces are those which 
are flat, whereas curved manifolds with (+ + ++) signature are known as Riemannian spaces. However, the 
terminology "Euclidean" action which has been taken over from flat space quantum field theory seems to 
have stuck, despite the fact that we are of course no longer dealing with K 4 . 

The action for fermi fields in ordinary quantum field theory is not positive-definite, but that is not a 
problem since one can treat them as anticommuting quantities so that the path integral over them converges. 
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wavefunction \1/ of the universe on a hypersurface, E, with intrinsic 3-metric, hij, 
and matter configuration, <£>, to be defined 38 ' 39 by the functional integral 



where the sum is over a class of 4- metrics, g^, and matter configurations, cf>, which 
take values hij and $ on the boundary E. Alternative definitions of the wavefunc- 
tion have been proposed. In particular, Linde 69 has argued that one should Wick 
rotate with the opposite sign, i.e., t — > +ir instead of t — > —it as above, lead- 
ing to a factor e +/ instead of e _/ in (3.11). Alternatively, one could stick with a 
Lorentzian path integral 70 , with e lS instead of e -7 in (3.11). In any case, in or- 
der to achieve convergence of the path integral® it is necessary to include complex 
manifolds in the sum, which somewhat obscures the distinctions between these al- 
ternative proposed definitions of The real distinction between the alternative 
definitions arises when it comes to imposing boundary conditions, thereby restrict- 
ing the 4-manifolds included in the sum in (3.11). For example, one could view the 
Euclidean sector as being the appropriate sector of the quantum theory in which an 
"initial" boundary condition on \1/ should be imposed, which would make (3.11) the 
natural starting point, as is the case for the no-boundary proposal 38 ' 39 . Alternative 
boundary conditions would favour the Lorentzian path integral 70 . 

The path integral definition of the wavefunction (3.11) is consistent with the 
earlier definition based on canonical quantisation to the extent that wavefunctions 
defined according to (3.11) can be shown to satisfy the Wheeler-DeWitt equation 
(3.8) provided that the action, the measure and the class of paths summed over are 
invariant under diffeomorphisms 71 . 

In the canonical quantisation formalism any particular solution to the Wheeler- 
DeWitt equation will depend upon the specification of boundary conditions on 
the wavefunction. In the path integral formulation the particular solution of the 
Wheeler-DeWitt equation that one obtains will similarly depend on the contour of 
integration chosen in superspace, and the class of 4-metrics one sums over in (3.11). 
Unfortunately it is not known how the choice of contour and class of paths prescribes 
the boundary conditions on the wavefunction in the general case, although it can 
be answered for specific models. The question of boundary conditions is naturally 
of prime importance for cosmology, and we shall return to this question in §4. 

3.4- Minisuperspace 

In practice to work with the infinite dimensions of the full superspace is not 
possible, at least with the techniques that have been developed to date. One useful 

® Linde's suggested modification [69] to (3.11) yields a convergent path integral for the scale factor, 
which is all that one needs in the simplest minisuperspace models, but does not lead to convergence if one 
includes matter or inhomogeneous modes of the metric. 




M 



(3.11) 
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approximation therefore is to truncate the infinite degrees of freedom to a finite 
number, thereby obtaining some particular minisuperspace model. An easy way to 
achieve this is by considering homogeneous metrics, since as was observed earlier 
for each point xGE there are a finite number of degrees of freedom in superspace. 

The results we shall obtain by this approach will be somewhat satisfying in that 
they do appear to have some predictive power. However, the truncation to minisu- 
perspace has not as yet been made part of a rigorous approximation scheme to full 
superspace quantum cosmology. As they are currently formulated minisuperspace 
models should therefore be viewed as toy models, which we nonetheless hope will 
capture some of the essence of quantum cosmology. Since we are simultaneously 
setting most of the field modes and their conjugate momenta to zero, which violates 
the uncertainty principle, this approach might seem rather ad hoc. However, in clas- 
sical cosmology homogeneity and isotropy are important simplifying assumptions 
which do have a sound observational basis. Therefore it is not entirely unreason- 
able to hope that a consistent truncation to particular minisuperspace models with 
particular symmetries might be found in future*. 

Let us thus consider homogeneous cosmologies for simplicity. Instead of having 
a separate Wheeler-DeWitt equation for each point of the spatial hypersurface, E, 
we then simply have a single Wheeler-DeWitt equation for all of S. The standard 
FRW metrics, with 



are of course one special example. In that case our model with a single scalar field 
simply has two minisuperspace coordinates, {a, $}, the cosmic scale factor and the 
scalar field. Many more general homogeneous (but anisotropic) models can also 
be considered. Indeed all such models can be classified^ and apart from the FRW 
models the other cases of interest are: (i) the Kantow ski- Sachs models 7 ' ' ,78 , which 
have a 3-metric 



and thus three minisuperspace coordinates, {a, 6, $}; and (ii) the Bianchi models. 

The Bianchi models are the most general homogeneous cosmologies with a 3- 
dimensional group of isometries. These groups are in a one-to-one correspondence 
with 3-dimensional Lie algebras, which were classified long ago by Bianchi 79 . There 
are nine distinct 3-dimensional Lie algebras, and consequently nine types of Bianchi 
cosmology. The 3-metric of each of these models can be written in the form 



h i jdx i dx j = a 2 (t) 



dr 2 



+ r 2 (dO 2 + sin 2 Odcp 2 ) , k = -1,0,1, (3.12) 



1 — kr 2 



/lijdxW = a 2 (t)dr 2 + b 2 (t) (d6 2 + sin 2 9dip 2 ) , 



(3.13) 



hijdx l dx J = hij{t)uj l <g> uj j , 



(3.14) 



* 



For some discussions of the validity of the minisuperspace approximation see [72-75]. 
See [76] for a review. 



18 



where u> 1 are the invariant 1-forms associated associated with the isometry group. 
The simplest example is Bianchi I, which corresponds to the Lie Group 1R 3 . In that 
case we can choose uj 1 = dx, uJ 2 = dy, and a; 3 = dz, so that 

h ij dx i dx j = a 2 (t)dx 2 + b 2 (t)dy 2 + c 2 (t)dz 2 , (3.15) 

and the minisuperspace coordinates are {a, 6, c, $}. If we take the spatial directions 
to be compact such a universe will have toroidal topology. In the special case that 
a(t) = b(t) = c(t) we retrieve the spatially flat (k = 0) FRW universe. 

The most complicated, and possibly the most interesting, Bianchi model is 
Bianchi IX, associated to the Lie group SO (3, M). In this case the invariant 1-forms 
may be written as 

uj 1 = — sin ipd9 + sin 9 cos ip dip, 

uj 2 = cos ip d6> + sin 9 sin if; dip, (3.16) 
u> 3 = cos 9 dip + dip, 

in terms of the Euler angles, (ip, 9, ip), on the 3-sphere, S 3 . The spatial sections of 
the geometry resulting from (3.14), (3.16) have the topology of S 3 , but are only 
spherically symmetric in the special case that /i 11 (t) = h 12 (t) = .../i 33 (t), which 
corresponds to the k = +1 FRW universe. Geometrically the spatial hypersurfaces 
can thus be thought of as squashed, twisted 3-spheres [see Fig. 3]. Bianchi IX has 
played an important role in classical cosmological studies of the initial singularity - 
it is the basis of the so-called "mixmaster universe" 80 ' 81 . As a classical dynamical 
system Bianchi IX is extremely interesting because it appears to be chaotic, but 
only just on the verge of being so. Over the years there has been much debate as 
to whether Bianchi IX is or is not chaotic, and this seems to have been recently 
resolved by an explicit demonstration that it is not integrable -82 . The correspond- 
ing minisuperspace model will have six independent coordinates in addition to the 
scalar field coordinate. 




(S* -^squash^- twist 



Fig. 3: Schematic geometry of spatial hypersurfaces in the Bianchi IX universe. 

' Technically speaking, what has been proved is the failure of integrability in the Painleve sense. While 
this does not guarantee the existence of chaotic regimes, it does make their existence extremely probable. 
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Let us now consider the minisuperspace corresponding to an arbitrary homoge- 
neous cosmology. We will assume that the minisuperspace is of dimension n, and 
will denote the minisuperspace coordinates by {q A }. Since Af l = by assumption, 
it follows from the definitions (2.5) and (3.4) that 



Q iiH KM = 4v^AT 2 (KijK^ - K 2 ) , 
and consequently the Lorentzian action (2.12) now takes the form 



S= dt 



2Af 



where 



Q AB dq A dq B = J d 3 x ^g^ kl 5h t3 5h kl + Vh5<S>5<S> 
is the minisupermetric, which is now of finite dimension, n, and 
U = J d 3 xVh ^ (~ 3r + 2A ) + V ( $ ) 



(3.17) 

(3.18) 
(3.19) 



(3.20) 



The action (3.18) is simply equivalent to that for a "point particle" moving in a 
"potential" U. Variation of (3.18) with respect to q A thus yields a geodesic equation 
with force term, 



l_d_ 

JfJt 



q_ 

Af 



+ ^r A 

M 2 i:: 



B .C n J 

q q = -G 



: dU 

dq B 



(3.21) 



where V A gc are the Christoffel symbols determined from the minisupermetric, while 
variation of (3.18) with respect to M yields the Hamiltonian constraint 



^ AB m Q +U(q) = 0. 



(3.22) 



The general solution to (3.21), (3.22) will have 2n — 1 independent parameters, 
one of which is always trivial in the sense that it corresponds to a choice of origin 
of the time parameter. In studying any particular minisuperspace model we must 
take care to check that what we have done above is consistent, as it does not always 
follow that substituting a particular ansatz into an action before varying it will yield 
the same result as substituting the same ansatz into the field equations obtained 
from variation of the original action. Eqs. (3.21) and (3.22) should correspond 
respectively to the (ij) and (00) components of the original Einstein equations, 
while the (Oi) equation is trivially satisfied in the present case. 

Quantisation is greatly simplified because now that our configuration space 
is finite-dimensional we are effectively dealing with the quantum mechanics of a 
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constrained system. The canonical momenta and Hamiltonian are respectively given 
by 



71, 



dL 



0q A M 



H = n A q A -L=M lg AB n A n B +U(q) 



J\fH. 



(3.23) 
(3.24) 



The 7t are related to the canonical momenta (2.13)-(2.16) defined earlier by inte- 
gration over the 3-volume of the hypersurfaces of homogeneity in (3.19). In terms 
of the new variables the action (3.18) and Hamiltonian constraint (3.22) are respec- 
tively 



S 



dt 



71 aQ 



\G AB n A n B +U(q) 



= 0. 



Under canonical quantisation (3.26) yields the Wheeler-DeWitt equation 

H* = [-lV 2 +U(q)} = 0, 



where 



V = 



1 



^GG AB d r 



(3.25) 
(3.26) 

(3.27) 
(3.28) 



is the Laplacian operator of the minisupermetric. In arriving at (3.27) we have made 
an explicit "natural choice" of factor ordering 33 ' 52 in order to accommodate the 
factor-ordering ambiguity. This choice is favoured by independent minisuperspace 
calculations of the prefactor, using zeta function regularisation and a scale invariant 
measure, which can then be related to factor ordering dependent terms through a 
semiclassical expansion of the Wheeler-DeWitt equation 83 . 

An alternative "natural choice" 35 > 84 ~ 86 of factor ordering would yield a confor- 
mally invariant Wheeler-DeWitt equation, 



1^9 n — 2 

2 8(n 



* = 0, 



(3.29) 



where TZ is the scalar curvature obtained from the minisupermetric. 



3.5. The WKB approximation 

In view of the difficulties associated with solving the Wheeler-DeWitt equation 
in general, the best we can realistically hope for in many minisuperspace models is 
to look for appropriate approximate solutions in the semiclassical limit 3 ", in which 

*~J]* n = J^e -1 ", (3.30) 

n n 

^ The semiclassical limit of the full superspace Wheeler-DeWitt equation has been treated more formally 
by a number of authors. See, e.g., [87], [88] and references therein. 
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where the sum is over saddle points of the path integral, the A n being appropriate 
(possibly complex) prefactors. In general we might expect to find regions in which 
the wavefunction is exponential, \& ~ e -7 , and regions in which it is oscillatory, 
~ e lS . The latter could be viewed as the wavefunction of a universe in the 
classical "Lorentzian" or "oscillatory" region, while the former would correspond 
to a universe in a classically inaccessible "Euclidean" or "tunneling" region. As 
has already been mentioned, the sum in (3.30) will in general contain a number of 
saddle points with an action, J n , which is neither purely real nor purely imaginary. 

Our own universe is of course Lorentzian at late times, and therefore the only 
minisuper space models which can be of direct physical relevance are those for which 
the Wheeler-DeWitt equation does possess approximate solutions of the oscillatory 
type. Approximate solutions of this type can be obtained by performing a WKB 
expansion, for which purpose it is necessary to restore h in the minisuperspace 
Wheeler-DeWitt equation (3.27). If we assume that each component \l/ n satisfies 
(3.27) separately, then 

= H^ n = [-\h 2 V 2 +U] A n ^ lh 
= e-Wft | [_i ( Wn ) 2 + w] A n + h [VI n ■ VA n + \A n V 2 I n ] + 0(ft 2 )} , 

(3-31) 

where the dot implies contraction with the minisupermetric Q . The 0(h°) and 
0(H) terms give two equations for I n and A n . If we decompose I n into real and 
imaginary parts according to I n = R n — iS n then the real and imaginary parts of 
the 0(ft°) term in (3.31) give 

-\ (VR n ) 2 + \ ( VS n ) 2 + U = 0, (3.32) 
VR n • VS n = 0. (3.33) 

Provided that the imaginary part of the action varies much more rapidly than the 
real part, i.e., (VR n ) 2 (VS n ) 2 , then (3.32) is the Lorentzian Hamilton- Jacobi 
equation for S n : 

i n AB dS n dS 
dq A dq 

Comparison of (3.34) with (3.27) suggests a strong correlation between coordinates 
and momenta, and invites the identification 

- |5- (3.35) 

If we differentiate (3.34) w.r.t. q c we obtain 

i r AB dS n dS n | n AB dS n d 2 S n | dU _ 

29 ' c dq A dq B+y dq A dq B dq c + dq c ~ U - ^ J 



2^ SS+^)=0- (3.34) 
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If we define a minisuperspace vector field 



d _ AB dS n d 

dS = g WW (3 ' 3 ) 

then combining (3.35), (3.36) and (3.37) we obtain 

^ + + J£ = 0, (3-38) 

which after raising indices is the same geodesic equation (3.21) obtained earlier 
provided we identify the parameter, s, with the proper time on the geodesies. 

We can now solve the equation given by the 0(h) term of (3.31). Since | VR n | <C 
|VS n | it follows that the terms involving R n can be neglected, and thus 

VS n • VA n = G AB ^ 3 -^ = ^ = -\A n V 2 S n , (3.39) 
oq oq ds z 

which may be readily integrated. We thus obtain a first-order WKB wavefunction 

tf n = C n exp ^zS n - \ J ds V 2 S n ^) , (3.40) 

where C n is an arbitrary (complex) constant to be appropriately normalised, and 
we have reverted to natural units in which h = 1. 

The wavefunction (3.40) could be considered to be the analogue of the wave- 
function for coherent states in ordinary quantum mechanics, 

i|>„(x, *) = Cne^" x exp ( ~ ( *~^ n(t))2 ) , (3.41) 

which describes a wave packet which is "peaked" about a classical particle trajec- 
tory, x n (£), and which thus roughly "predicts" classical behaviour. This becomes 
problematic, however, if we consider a superposition, i|> = ^2 n i\> n , of such states 
since interference between different wave packets will in general destroy the classical 
behaviour. In order to interpret the total wavefunction as saying that the particle 
follows a roughly classical trajectory, x n (t), with probability |c n | 2 , it is necessary 
that a decoherence mechanism should exist which renders this quantum mechanical 
interference negligible 89 . 

The issue of quantum decoherence is clearly also of great importance in quantum 
cosmology, since in order to interpret \1/ in (3.30) in a similar fashion a similar mech- 
anism must exist. It has been recently argued, furthermore, that decoherence is a 
necessary feature of the WKB interpretation of quantum cosmology, since without 
decoherence the existence of chaotic cosmological solutions would lead to a break- 
down of the WKB approximation 90 . This is analogous to similar problems with 
the commutativity of the limits t — > oo and h — > in ordinary quantum mechanics 
when applied to chaotic systems. 

The issues involved in decoherence pose complex conceptual questions for the 
fundamentals of quantum mechanics itself, quite apart from the problems specific 
to quantum cosmology^ 1 . Here we will merely assume that such a mechanism exists, 

For further details see, e.g., [91]- [94] and references therein. 
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and we will take the view that \1/ can be considered to "predict" a classical spacetime 
if there exist WKB-type solutions (3.40), which yield a strong correlation between 
n A and q A according to (3.35). The sense in which the minisuperspace positions 
and momenta are "strongly correlated" can be made more precise through the use 
of quantum distribution functions, such as the Wigner function 91 ' 95 . By use of the 
Wigner function one may show 95 that wavefunctions of the oscillatory type, ^> ~ e lS , 
predict a strong correlation between coordinates and momenta, whereas wavefunc- 
tions of the type \1/ ~e _I , which are also typical minisuperspace solutions, do not. 
Such exponential wavefunctions can thus be considered as describing universes in 
a purely quantum "tunneling" regime, before the quantum to classical transition. 
We will interpret wavefunctions, \1/ ~e lS , as corresponding to classical spacetime, or 
rather a set of classical spacetimes as S is a first integral of the equations of motion. 

3.6. Probability measures 

Given a solution, to the Wheeler-DeWitt equation it is necessary to con- 
struct a probability measure in order to make predictions. One central question in 
quantum cosmology is how one should construct such a measure. 

The minisuperspace Wheeler-DeWitt equation (3.27) is a second-order equation 
very much like the Klein-Gordon equation in ordinary field theory, and it readily 
follows from (3.27) that the current 33 

a = -\i (^"V^ - tfVtf) (3.42) 

is conserved: V • 3 = 0. The similarity to the Klein-Gordon current extends to 
the fact that the natural inner product 33 constructed from 3 is not positive-definite 
and so gives rise to negative probabilities. In quantum field theory this is not a 
problem since one can split the wavefunction up into positive and negative fre- 
quency components which correspond to particles and anti-particles. However, as 
has already been mentioned there is no well-defined notion of positive frequencies in 
superspace on account of its lack of symmetries 96 . A further problem is that many 
natural wavefunctions would have zero norm with this definition. For example, the 
no-boundary wavefunction is real and gives 3 = 0- 

The similarity of ^ to the Klein-Gordon field has suggested to many people that 
one should turn ^ into an operator, ^, thereby introducing quantum field theory on 
superspace, or "third quantisation" . One then arrives at operators which create and 
annihilate universes. However, as we do not perform measurements over a statistical 
ensemble of universes it is not clear how we can arrive at sensible probabilities using 
such a formulation. 

The difficulties with the Klein-Gordon current of course led Dirac to introduce 
the Dirac equation, and it is worth mentioning that a similar resolution of the 
problem is available in supersymmetric quantum cosmology. In particular, one can 
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go to a theory which includes fermionic variables by considering quantum cosmol- 
ogy based on supergravity"* rather than the purely bosonic Einstein theory. The 
constraints of supergravity, which may be viewed as the Dirac square root of the 
constraints of general relativity 98 ' 99 , are reducible to first-order equations. Fur- 
thermore, this also translates into simplifications in homogeneous minisuperspace 
models - the appropriate constraint equation which determines the quantum evo- 
lution of the wavefunction can be considered to be the Dirac square root of the 
Wheeler-DeWitt equation 100-103 . As a result it is possible to construct 100 a Dirac- 
type probability density which is conserved by the equation Q^> = 0, where Q is 
the supercharge. 

Another alternative to the question of the probability measure is to use |*| 2 
directly as a probability measure 52 ' 104 , by defining the probability of the universe 
being in a region, A, of superspace by 



where *1 is the volume-element on superspace, * being the Hodge dual in the 
supermetric. This definition of a necessarily positive-definite probability density 
works very well for homogeneous minisuperspaces, for which the volume form * 1 is 
independent of x G E. This is perhaps not surprising since as was observed in §3.4 
the assumption of homogeneity reduces the problem to one of quantum mechanics, 
and |\1/| 2 is of course the probability density in conventional quantum mechanics. 

Problems with the definition (3.43) do arise since even in some simple examples 
the wavefunction is not normalisable, but instead (^1^) = oo. One further problem 
is that whereas in ordinary quantum mechanics |\1/| 2 describes a probability density 
in configuration space - i.e., the space of particle positions - in quantum cosmology 
the configuration space is (mini) superspace and time is implicitly contained in the 
(mini) superspace coordinates. These coordinates cannot therefore be thought of as 
the mere analogues of particle positions. As a result the recovery of the conservation 
of probability and the standard interpretation of the quantum mechanics for small 
subsystems is not necessarily straightforward in the approach based on (3.43), and 
may involve understanding some subtle questions about the role of time in quantum 
gravity. 

Ultimately a formulation such as (3.43) which is based on absolute probabili- 
ties may not be required since it is impossible to measure statistical ensembles of 
universes and thus all we can really test are conditional probabilities rather than 
absolute probabilities. For example, a relevant testable probability might be the 
probability, V(A\H), of finding \1/ in a region A of superspace given that \1/ started 
in another region £ of superspace. Page 2 ' 105 has explored the construction of condi- 
tional probabilities in quantum cosmology without the use of absolute probabilities. 

^ In fact, a naive first-order Hamiltonian formulation for the minisuperspaces of homogeneous cosmolo- 
gies was found early on [97], but until the development of supergravity there was no natural interpretation 
of the Dirac-type constraint equation obtained. 




(3.43) 
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Clearly the issues surrounding the choice of probability measure involve some 
deep conceptual problems which may perhaps get to the heart of the broader con- 
ceptual basis of quantum gravity. Such issues have been discussed by a number 
of authors in the context of quantum cosmology 2 ' 43 ' 104-107 and I will not address 
them here in any detail. 

For the purposes of examining how we might hope to make predictions from 
the proposed boundary conditions of §4.1,4.2 we shall merely consider quantum 
cosmology in the WKB limit. In this limit it follows from (3.40) that each of the 
components ^ n in (3.30) has a conserved Klein-Gordon-type current (3.42) given 
by 

a n ~ |-4 n | 2 VS n , (3.44) 

which flows very nearly along the direction of the classical trajectories. The current 
conservation law V A 3 n A = implies that 

dV = 3 n A dL A , (3.45) 

is a conserved probability measure on the set of trajectories with tangent VS n , 
where dL A is the element of a hypersurface, L, in minisuperspace which cuts across 
the flow and intersects each curve in the congruence once and only once. One 
finds that for a pencil of trajectories near the classical trajectory the probability 
density (3.45) is positive-definite. Vilenkin has argued 104 that positive-definiteness 
of the probability measure is really only required in the semiclassical limit, as this 
is the only limit in which we obtain a universe accessible to observation where the 
conventional laws of physics apply. Therefore given that the current (3.42) works 
in the WKB limit, this is all that is needed if we are content that "probability" 
and "unitarity" are only approximate concepts in quantum gravity. One can also 
show 2 ' 52 that the manifestly positive definition (3.43) yields essentially the same 
result as (3.45) in the WKB limit. 



3. 7. Minisuperspace for the Friedmann universe with massive scalar field 

Let us now apply our results to the particularly simple case of a homogeneous, 
isotropic universe with a single scalar field, with a potential which allows for infla- 
tionary behaviour. A quadratic potential is possibly the simplest example with this 
property, and thus has been much studied in quantum cosmology. 

For convenience we will introduce a numerical normalisation factor a 2 = 
k 2 /(3V) into the metric, where V is the 3-volume of the unit hypersurface - e.g., 
V = 2-7T 2 for the 3-sphere, k = +1. We are considering closed universes only, which 
requires making topological identifications for the k = 0, — 1 cases, so that V remains 
finite. In place of (2.1), (2.2) and (3.12) we then have a metric 



ds' 



a 2 \-M 2 dt 2 + a 2 (t) 



dr 2 
1 — kr 2 



+ r 2 (d9 2 + sin 2 6dp 2 ) 



(3.46) 
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We also take the scalar field action to be normalised by 

S mattCT = ^ / {—Q^d^ - ^jpj . (3.47) 

M 

Thus J\f — > aAf, a — > aa, $ = v/30/k and V = 3V/(2« 2 a 2 ) relative to our earlier 
definitions, and a, 4> and F are now dimensionless. As for the general homogeneous 
minisuper space discussed in the previous section, we use a gauge in which Mi = 0. 
From (3.46) it follows that 3 R = and Kij = —-^jj^hij, and thus the action 
takes the form (3.18) with a minisupermetric 

G AB dq A dq B = -ada 2 + a 3 d<p 2 , (3.48) 

and potential 

U=\ (a 3 V(<i)) - ka) . (3.49) 
Alternatively, it is sometimes useful to express (3.48) in terms of a conformal gauge 

Q AB dq A dq B = e 3a (-da 2 + d</> 2 ) , (3.50) 
= -(Auv)~ 1/A dudv, (3.51) 

where a = In a, or alternatively in null coordinates 



u =l e 2 ^ = \a\~ 2 \ 
v =Ie 2 ^ = \a 2 ^. 



(3.52) 



The canonical momenta are given by 



71 = U, 7t a = 7t^ = (3. 53) 

and the classical equations of motion yield 

,3^2 \ 

-ka + a 3 V J =0, (3.54) 

2ad> 2 

MdtKM)^-^^ (3 ' 55) 

Id ( 4>\ a4> ; 1 dV 
A7d7( v A7j +3 ^ + 2 d^ = ' (3 ' 56) 

which are equivalent to the geodesic equation (3.21) with metric (3.48) and potential 
(3.49). The lapse function is not of physical relevance classically since we can choose 
an alternate proper time parameter, dr = Afdt, or equivalently choose a gauge 
Af = 1 in (3.54)-(3.56) so that t is the proper time. One of these equations depends 
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on the other two by virtue of the Bianchi identity, as is always the case in general 
relativity. We thus effectively have two independent differential equations in two 
unknowns, one first order and one second order, or equivalently an autonomous 
system of three first order differential equations. The solution therefore depends on 
three free parameters, but as mentioned above one of these amounts to a choice of 
the origin of time which is of no physical importance. Thus there is a two-parameter 
family of physically distinct solutions. 

The classical solutions cannot be written in a simple closed form except for 
certain special values of k and the potential V ((/)), which unfortunately does not 
even include the quadratic potential^ V(<f>) = m 2 (f) 2 . The qualitative property of 
the solutions may nonetheless be determined by studying the 3-dimensional phase 
space. Instead of choosing a particular potential, however, let us suppose that we 
are in a region of the phase space for which V can be approximated by a constant. 
Such conditions are more or less met in the "slow-rolling approximation" 109 of in- 
flationary cosmology, in which <^ 6 and |V"/V| <C 9. In this approximation 
the dynamics is described by setting <fi ~ in (3.54) and (3.55), and setting <j> ~ 
in (3.56). In approximating the nearly flat region of the potential V by a cosmolog- 
ical constant we ignore the slow change of the scalar field determined from (3.56). 
We thus obtain a simplified model which possesses classical inflationary solutions, 
provided the constant V is chosen to be positive. This will serve as a useful test 
model for quantum cosmology. 

In the case that V is constant, eq. (3.56) integrates to give <p = Ca -3 , where 
C is an arbitrary constant, and it therefore follows that the Friedmann equation 
(3.54) can be written in terms of an elliptic integral in a 2 (r]), 

r, = \ r dZ = , (3.57) 

2 Jo Wz 3 - kz 2 + C 2 V ' 

where rj is the conformal time parameter defined by dry = a _1 dt. It is thus possible 
to express the general solution in terms of elliptic functions. Since the properties 
of such functions are not very transparent perhaps, we can alternatively plot the 
2-dimensional phase space - e.g., in terms of the variables <f> and a, as in Fig. 4 in 
order to understand the qualitative features of the solutions. There are four critical 

points, at (0, a) = | ^0, ±l/VV^j , (±l/s/2V, j, the first two being nodes, A±, 

which are endpoints for all values of the spatial curvature, k, and the latter two 
saddle points, B±, for the k = +1 solutions. 

The general solution for the spatially flat case (k = 0), which corresponds to 
the bold separatrices shown in Fig. 4, is given by 110 



a = C, 



1/3 

sinh | V Vt cosh ^VVt 



</>= | In 



tanhf 



2 



1/3 



(3.58) 



^ The Einstein equations for the FRW universe with a massive scalar field can be solved approximately 
[108] in various limits, however. 
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Fig. 4: The 2-dimensional phase plot of a — a/ a versus <j> for the simplified model 
with constant potential, V. 

in closed form, where C a and are arbitrary constants. At late times the solutions 
(3.58) have an exponential scale factor, a — > C a exp(Vv t) as t — > oo, and constant 
scalar field <p — > C^,. Furthermore, one can see from Fig. 4 that all the k = — 1 
solutions in the upper half-plane, and a number of the k = +1 solutions are also 
attracted to the point A + with a similar inflationary behaviour at late times. (The 
corresponding point A_ on the k = separatrix corresponds to the time-reversed 
solution, with an inflationary phase as t — > — oo.) The simplified universe corre- 
sponding to Fig. 4 of course is far from being the complete picture, as the model 
does not allow for any exit from inflation. However, Fig. 4 illustrates the typical 
situation that a given model will possess regimes with inflationary behaviour and 
regimes with non-inflationary behaviour. In Fig. 4 the k = +1 solutions to the right 
and left of the separatrices that pass through B± fall into the latter category, for 
example. The situation becomes even more involved when one considers the full 
3-dimensional phase space for some particular potential V ((/)). 

The case of the k = +1 solutions in Fig. 4 illustrates the general feature that 
classical dynamics are highly dependent on initial conditions. In order to obtain 
a sufficiently long inflationary epoch to overcome the problems mentioned in the 
Introduction, (of order 65 e- folds growth in the scale factor), the initial values of <fi 
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and (f> must be restricted to a particular region of the phase space. In particular, 
4> must be small initially. Classically, there is no a priori reason for one choice of 
initial conditions over any other choice, unless further ingredients are added. The 
degree to which inflationary initial conditions are preferred relative to other initial 
conditions - i.e., how probable is inflation? - is precisely the sort of question that 
we might therefore hope quantum cosmology could answer. 

It is possible to attempt to solve this question without resorting to quantum cos- 
mology. To do this one must construct a measure on the set of all universes 111 ' 112 , 
and then compare the number of inflationary solutions with a sufficiently long ex- 
ponential phase to the number of other solutions. Preliminary results 112 seemed to 
indicate that almost all models with a massive scalar field undergo a period of infla- 
tion. However, a more careful analysis 108 revealed that the answer is ambiguous, as 
both the set of inflationary solutions and the set of non-inflationary solutions have 
infinite measure. 

Alternatively, if as we expect the universe began in some sort of tunneling 
process or similar transition from a quantum regime, then we could expect the 
"initial" classical parameters to be determined, at least in a probabilistic fashion, 
from more fundamental quantum processes. The question of the most probable 
state of the universe is then pushed back a level and becomes: "what is a typical 
wavefunction for the universe?" 

In the context of the present minisuperspace model, therefore, we can proceed 
by quantising the Wheeler-DeWitt equation (3.54), to obtain 



= [-±V 2 +U] * 



a; 



d d 



d 2 



Q da Q da d<j) 2 



^ - ka + a 3 V((f)) 



1 -3a 
2 e 



O 2 



d 2 



- ke 4a + e ba V((f>) 



6ai 



da 2 dcj) 2 



(3.59) 



(4m;) 1/4 




F) 2 h 



in terms of the various sets of coordinates given earlier. In general boundary con- 
ditions will have to specified in order to solve (3.59). However, we can consider the 
approximate form of the WKB solutions without considering boundary conditions 
for the time being. 

We will confine ourselves to regions in which the potential can be approximated 
by a cosmological constant, as in the analysis of Fig. 4, so that we can drop the 
term involving derivatives with respect to in (3.59), thereby obtaining a simple 
1-dimensional problem which is amenable to a standard WKB analysis. The first 
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order WKB wavefunction (3.40) which solves (3.59) in this approximation is 



#(a,0) 



3(0) 



a(a 2 V(<P)-k) 1/4 



k a(k- a 2 V((j))) 



1/4 



exp 



exp 



±i 



3V(0) 
±1 

WW) 



(a 2 V(<P)-k) 3/2 
(fc-aV(0)) 3/2 



a 2 V > k, (3.60) 



a 2 V < k. (3.61) 



If V is positive, as was assumed above, then oscillatory type solutions will thus exist 
for large values of the scale factor, while the exponential type solutions will exist 
for small values of the scale factor if k = +1. 

The oscillatory solutions are of the form \l/~e lS (neglecting the prefactor), 
where S satisfies the Hamilton- Jacobi equation (3.34). Comparing this to the Hamil- 
tonian constraint (3.26) we find a strong correlation (3.35) between momenta and 
coordinates. For large scale factors, a 2 V ^> \k\, so that S ~ ±^a 3 VV. In this limit 
(3.23), (3.35) and (3.50) thus yield 



OS 

oa 

ds 

n * = d} 



a 



~ o, 



(3.62) 



which correspond in fact to the inflationary points, A±, of Fig. 4. The oscillatory 
wavefunction thus "picks out" classical inflationary universes. 

Since the minisupermetric (3.50) is conformal to 2-dimensional Minkowski space 
in the coordinates (a, 0), it is convenient to represent it by a Carter-Penrose con- 
formal diagram (see Fig. 5). In each case we plot (p — q) horizontally and (p + q) 
vertically, where tanp = a + 4>, and tang = a — 4>. The boundary consists of 
points corresponding to past timelike infinity, x~ = {(a, 0) | a = 0, finite }, fu- 
ture timelike infinity, i + = {(a,</>) |a = oo,0 finite }, left and right spacelike in- 
finity, i° L R = {(a, 0) | a = finite , = ±oo}; and past and future null boundaries, 
3~ = {(a, 0) | a = 0, = ±00} and J+ = {(a, 0) | a = 0, = ±00}. In each case 

the subscript L (left) is associated with — > —00, and the subscript R (right) 
with — > +00. The approximate region for which oscillatory WKB solutions exist 
is shown in Fig. 5(a,b) for the approximate minisuperspace with a cosmo logical 
constant, in Fig. 5(c,d) for V(<p) = m 2 (j) 2 , and in Fig. 5(e,f) for potentials, V(0), 
typically found in higher-derivative gravity theories and in string theory with su- 
persymmetry breaking. 

Naturally it is of interest to know whether the inflationary WKB wavefunctions 
are typical solutions to the Wheeler-DeWitt equation. To determine a typical wave- 
function for the universe, we need to make a choice of boundary conditions for \1/ 
when solving (3.59). 
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Fig. 5: Conformal diagrams of the 2-dimcnsional minisupcrspace. The region where 
oscillatory WKB solutions exist, as given by the rough criterion a 2 V > 1, is shaded 
for various potentials: (a) V = 0.25 (const); (b) V = 4 (const); (c) V = 0.25</> 2 ; 

(d) V = 25(f) 2 ; (e) V = [l - e"^] 2 with / = 1.5; (f) V = 4sinh 2 0exp [-/e- 2 *] 
with / = 0.1. 
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4. Boundary Conditions 

The specification of boundary conditions for the Wheeler-DeWitt equation may 
seem a disappointment, as it might appear that we are just replacing an arbitrary 
initial choice of parameters which describe the classical evolution of the universe 
by an arbitrary initial choice of parameters which describe its quantum evolution. 
However, if quantum mechanics is a universal theory then it must have applied at 
the earliest epochs of the existence of the universe, in which case it is natural that 
the quantum dynamics precedes the classical dynamics. This justifies a quantum 
boundary condition for the universe as being more fundamental than a classical one. 
In any case, the only alternative to choosing quantum boundary conditions would 
be that mathematical consistency might be enough to guarantee a unique solution 
to the Wheeler-DeWitt equation, as DeWitt originally hoped 33 . If the experience 
gained from the study of minisuperspace models translates to superspace, then this 
would not appear to be the case, however. 

The question naturally arises as to whether there should be some natural bound- 
ary condition, which once and for all determines the quantum evolution of the uni- 
verse at early times, or alternatively whether the nature of the quantum dynamics 
might be somewhat indifferent to such choices. Deep conceptual problems are in- 
volved in trying to make headway with this question. Unlike other situations in 
quantum physics, where boundary conditions are readily specified by the symme- 
try of particular problems, such as spherical symmetry in the case of the hydrogen 
atom, the origin of the universe poses a situation in which all intuition must be 
abandoned and we can at best proceed on aesthetic grounds alone. 

Having made a choice of boundary condition, we can of course solve the Wheeler- 
DeWitt equation and study the physical consequences for the evolution of the uni- 
verse. However, without some additional principle we should by rights study many 
different boundary conditions before we can begin to have any confidence about the 
predictions made. To arrive at a principle which would circumvent this problem is 
an immense challenge: it would more or less amount to an additional law of physics 
which must be appended to the others which describe the quantum evolution of the 
universe. The situation might be considered to be the same as trying to describe 
the phase transition from gas to liquid if all the physical phenomena that we knew 
about related to the gaseous phase only. The universe appears to have undergone 
a phase transition when it was formed, but the only experience we have available 
involves the "after" state of the universe alone. 

Progress can of course only be made by attempting to define natural boundary 
conditions for the wavefunction of the universe, and examining the consequences. 
This became an important activity in the 1980s. I will only discuss the two most 
studied boundary condition proposals, the "no-boundary proposal" and the "tun- 
neling proposal". However, other proposals have been put forward, including the 
"all possible boundaries proposal" of Suen and Young 113 and the "symmetric initial 
condition" of Conradi and Zeh 114 ' 115 . 



33 



4-1. The no-boundary proposal 

The proposal of Hartle and Hawking 38 ' 39 is that one should restrict the sum in 
the definition of the wavefunction of the universe (3.11) to include only compact 
Euclidean 4-manifolds, M, for which the spatial hypersurface E on which ^> is 
defined forms the only boundary, and only matter configurations which are regular 
on these geometries. The universe then has no singular boundary to the past, as is 
the case for the standard FRW cosmology. The sum (3.11) thus includes manifolds 
such as those shown in Fig. 6, but not those shown in Fig. 2. As Hawking 39 puts 
it: the boundary conditions of the universe are that it has no boundary. 




Fig. 6: Geometries allowed by the Hartle-Hawking no-boundary proposal. 



Intuitively, what Hartle and Hawking had in mind in formulating this proposal 
was to get rid of the initial singularity by "smoothing the geometry of the universe off 
in imaginary time" . For example, whereas a surface with \fh = would be singular 
in a Lorentzian signature metric, this is not necessarily the case if the metric is 
of Euclidean signature, as can be seen from the example of S 4 shown in Fig. 7. 
Ideally, the no-boundary proposal should tell us what initial conditions to set when 
we take manifolds with \fh — > or any similar limit consistent with the proposal. 
If the limit is taken at an initial time r = 0, the no-boundary proposal would 
lead to conditions on hij(x,0), <f>(x,Q) and their derivatives. In practice, quantum 
cosmology is rarely studied beyond the semiclassical approximation, in which ^> ~ 
Ae~ Icl , where I c \ is the classical (possibly complex) action evaluated along the 
solution to the Euclidean field equations. In the semiclassical approximation one 
therefore works only with boundary conditions on the metric and matter fields 
which correspond to the no-boundary proposal at the classical level. In particular, 
we demand: (i) that the 4-geometry is closed; and (ii) that the saddle points of the 
functional integral correspond to regular solutions of the classical field equations 
which match the prescribed initial data on E. 

One question which is not explained by the no-boundary proposal is the choice 
of a contour of integration for the path integral. As was mentioned earlier, the path 
integral over real Euclidean metrics does not converge, and thus it is necessary to 
include complex metrics to make the path integral converge. Such metrics will gen- 
erally include ones which are neither truly Euclidean nor truly Lorentzian, and thus 



34 




Fig. 7: Slicing a 4-sphere of radius R embedded in flat 5-dimensional space: 

(a) a surface \x 5 \ — a < R intersects the 4-sphere in a 3-sphere of non-zero radius; 

(b) when \x 5 \ —R the 3-sphere shrinks to zero radius but there is no singularity of 
the 4-geometry. 



the naive picture of a compact Euclidean geometry sewn onto a Lorentzian one, 
which is suggested by "smoothing the geometry of the universe off in imaginary 
time", is not completely accurate. In general, one might expect a truly complex 
metric to interpolate the Euclidean and Lorentzian ones, and in general the ini- 
tial geometry might be only approximately Euclidean and the final geometry only 
approximately Lorentzian 116 . Unfortunately the criteria for achieving convergence 
of the path integral do not single out a unique contour of integration, and the no- 
boundary proposal does not appear to offer any further clues as to how the contour 
should be chosen 117 . 




Fig. 8: Euclidean solutions which correspond to matching a given 3-sphere hy- 
persurface, S, to a 4-sphere which is: (a) less than half filled; (b) more than half 
filled. 

Non-uniqueness of the contour of integration is already a problem in the simplest 
conceivable non-trivial minisuperspace model, namely the k = +1 FRW universe 
with a cosmological term and no other matter - the "de Sitter minisuperspace" 
model. At the semiclassical level one can calculate e~ 7cl by the steepest descents 
method. Hartle and Hawking discussed this in their original "no boundary" paper 38 , 
and argued heuristically that one particular saddle point would yield the dominant 



35 



contribution to the path integral, namely the saddle point corresponding to the 
classical Euclidean solution which matches the S s hypersurface E to a less than 
half filled 4-sphere, rather than the solution which matches S 3 to a more than 
half filled 4-sphere (see Fig. 8). However, their argument did not stand up to a 
more rigorous analysis. Halliwell and Louko 66 found a means of evaluating the 
minisuperspace path integral exactly, and thereby explicitly determined convergent 
contours. They showed that this simple model possesses inequivalent contours for 
which the path integral converges. These pass through different saddle points and 
lead to different semiclassical wavefunctions, There are thus many different 
no-boundary wavefunctions, each corresponding to a different choice of contour. 
The problem persists in more complicated models 67 ' 68 . Since different no-boundary 
wavefunctions could lead to different physical predictions, the ambiguity associated 
with the choice of contour would appear to be the most significant problem with 
the no-boundary proposal which still needs to be resolved. 

Let us return to the minisuperspace example of §3.7 and examine the implica- 
tions of imposing the Hartle-Hawking boundary condition. For simplicity we will 
specialise to k = +1 models. The other values of k have also been discussed in the 
literature. 2 ' 105 ' 118 

The minisuperspace no-boundary wavefunction for the k = +1 models is given 

by 

# HH [a,0] = J ' DaSct^A/'e-^'^ 1 , (4.1) 



where 1 



™T f 

7=1 / drAT 



-a / da\ 2 a 3 / d<±A 2 3ta 
- + ttk I -r~ -a+a 3 V 



(4.2) 



Af 2 \drj Af 2 \dr 
The integral is taken over a class of paths which match the values 

a(r f ) = a, <J>(r/) = 0, (4.3) 

on the final surface, and the origin of the Euclidean time coordinate, r, has been 
chosen to be zero. 

There are two approaches we can take to solving the Wheeler-DeWitt equa- 
tion (3.59): either (i) attempt to interpret the Hartle-Hawking boundary condition 
directly in terms of boundary conditions of ^ on minisuperspace; or (ii) take a 
saddle-point approximation to the path integral (4.1). For consistency these two 
approaches should agree. 

Let us first consider (3.59) directly. Hawking and Page 52 ' 119 have argued that 
one can approximate the Hartle-Hawking boundary condition in minisuperspace by 
saying that in an appropriate measure one should have \& = 1 when a — * with (f> 



A I shall use a different font to distinguish the minisuperspace coordinates in the functional integral, 
(a, 4)), from their boundary values, (a, 0), in the 3-geometry on the hypersurface E. 
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regular, and ^ = 1 also along the past null boundaries, in order to provide sufficient 
Cauchy data to solve (3.59) everywhere in the (a, 0) plane. It is possible to solve 
(3.59) exactly for a massless scalar field 86 ' 120 , (i.e., V = 0), but exact solutions 
are not known for the massive scalar field. Approximate solutions can be found in 
various regimes 52 ' 119 ' 121 , however. 

Firstly, since ^ must be regular as a — > 0, we see that \1/ must be independent of 
4> in this limit, ^| ~ 0, in order to overcome the divergence of the a~ 3 factor (using 
coordinates (a, (p)) in (3.59). Furthermore, for a 2 V <C 1 the approximation V ~ 
is a good one M , and for k = ±1 (3.59) becomes a Bessel equation in the variable 
\c? ■ For k = +1, which is the case of interest to us here, we therefore obtain the 
solution 52 ' 119 

^I (±a 2 ), (4.4) 

where l (z) = J2^Lo j^ryz (f ) 2 ™ is the zero order modified Bessel function, and the 
normalisation has been fixed to satisfy the boundary condition f -> 1 as o -> 0. 
For large a (with a 2 V <ti 1), 

*~-^exp(±a 2 ) [l + 0(a" 2 )], (4.5) 

\/ 7T Q. 



i.e., the wavefunction is of exponential type. It therefore agrees with the WKB 
approximation (3.61) in the limit a 2 V <^ 1 for large a provided we take the (— ) 
solution of (3.61) with normalisation 

c «^^(4)' (46) 

Let us now consider the limit V((f>) ^> 1 but avoid regimes in which the 4> dependence 
of (3.59) is significant by assuming that V is approximately constant, as in the 
approximation of Fig. 4, so that the <p derivative can still be neglected. The term 
involving the spatial curvature k in (3.59) is now negligible compared to the last 
term and can also be neglected, yielding the solution 52 ' 119 

¥~c(0)J o (|aVn (4-7) 

where J (^) = Zl^Lo T^Tp (f) 2 " 1S the zero order ordinary Bessel function. Since 
the spatial curvature term in (3.59) is always dominant for a — > 0, the approximation 
which led to (4.7) no longer applies in that limit, so the factor c cannot be normalised 
by the boundary condition \l/(0, 0) = 1. For large a, however, 

*~7fer os ( s -i)' (48) 



1X1 We assume that V(cj>) grows less strongly than e 6 '* 1 ' as \(f>\ — > oo, i.e., | V'/V| < 6, so that this 
approximation remains valid for arbitrarily large \(f>\. 
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where S = ^a 3 W, which is a superposition of the two oscillatory WKB modes 

(3.60) for large S. Using the WKB connection formula to match the (— ) solution of 

(3.61) to the oscillatory region, we find agreement with the asymptotic limit (4.8), 
provided 

^2tT /T (I 




c=^ T C = V|exp^j. (4.9) 

The ^dependent corrections to these wavefunctions, which result from perturba- 
tions in directions in which 0^0 have been discussed by Page 2 ' 119 ' 122 . 

Let us now consider the alternative method of determining \& by making a 
saddle-point approximation to the path integral. In the semiclassical approximation 
the wavefunction takes the form 

tf~exp(-J cl (a,0)), (4.10) 

where I cl denotes the Euclidean action (4.2) evaluated at a classical (possibly com- 
plex) solution to the Euclidean field equations 

(£) , - i+a * v( «- - (4 - n) 

11 J_da d* _ ,dV = 

MAt\MAt) Afadrdr 2 d<|> ' K ' 

which follow from (3.54)-(3.56) by replacing t — > —it. We will henceforth restrict 
ourselves to the gauge in which ^ = 0, i.e., Af = const, in which case the functional 
integral over Af in (4.1) must be replaced by an ordinary integral. 

The Hartle-Hawking boundary condition demands that 



= 0. (4.14) 

o 



To see this consider the Euclidean 4-metric which is given by 

ds 2 = AT 2 dr 2 + a 2 (r)d0 3 2 . (4.15) 

The Hartle-Hawking boundary condition requires that we close this 4-geometry in a 
regular fashion as r — > 0. This is achieved if a~ ±A/Y as r — > 0, since (4.15) is then 
the same as the metric of the 4-sphere in spherical polar coordinates in this limit. 
This suggests that we demand a(0) = and jf = ±1- However, the second 
condition is guaranteed by the constraint equation (4.11) if the first condition is 
imposed. The second condition of (4.14) is obtained by noting that (4.13) will give 



a regular solution for 4> in the limit r — > only if ^ 
of (4.13) diverges otherwise. 



0, since the middle term 



38 



If as before we make the simplifying assumption that V can be approximated 
by a cosmo logical constant, and -£p ~ 0, then it follows that there exist solutions 
satisfying the boundary conditions (4.3)and (4.14), which may be written 

sin(.A/VV>/) 

where 

sm 2 (tf VVTf) = a 2 V, (4.17) 

which follows from solving the constraint (4.11). If a 2 U < 1 then (4.17) will give 
an infinite number of real solutions for the constant A/", which may be conveniently 
parametrised 

J\f = Af n ± = —fL — \{n + |)7r ± cos~ 1 (av / V F )j , n G Z, (4.18) 
\JV Tf L J 

with cos _1 (aV / V F ) taken in the principal range (0, ^). Substitution of (4.18) in (4.16) 
then gives 

a(r) ~ (-l) n y- 1/2 sin(Arv / Fr). (4.19) 

In addition to the real solutions (4.18) of (4.17), there will also be complex solutions 
if a 2 V > 1. 

Using the solution (4.19) it is straightforward to evaluate the classical action 
(4.2). For n = 0, for example, we find 



I± = 



■1 r. . /- 2 T ^/.\\3/2' 



3^(0) 



lifl-oV^))' . (4.20) 



If we substitute (4.19) into (4.15) we obtain the metric of the 4-sphere, and thus the 
classical solutions correspond to matching a given 3-sphere to 4-sphere(s). Further- 
more, for n = we find ^| = =fA/\/1 — a 2 V, so that the (— ) solution of (4.18) 

and (4.20) has ^| > 0, which corresponds to matching the 3-geometry to a less 

than half filled 4-sphere (Fig. 8(a)), while the (+) solution similarly corresponds to 
the more than half filled 4-sphere case (Fig. 8(b)). Values of n > would appear to 
give cases in which the 4-geometry pinches off to zero a number of times and then 
"bounces back" resulting in linear chains of contiguous 4-spheres 123 ' 124 . However, 
this interpretation of the saddle points as corresponding to universes with bounce 
solutions does not appear to remain valid 116 once one considers complex solutions 
to the field equations (4.11)-(4.13). 

As was discussed earlier the no-boundary condition does not prescribe a unique 
contour of integration, and thus it is not completely clear which of the above sad- 
dle points should be included in the semiclassical wavefunction (3.30). Halliwell 
and Hartle have shown 117 that the points with n < 0, which have a negative lapse 
function, lead to difficulties with the recovery of quantum field theory in curved 
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spacetime from quantum cosmology, and thus one might hope that these points 
should be avoided in the contour of integration. However, no clear grounds present 
themselves for omitting the other saddle points. For the purpose of making predic- 
tions from the no-boundary proposal we shall therefore make a choice by assuming 
that the contour is such that the solution corresponding to the less than half filled 
4-sphere, with action 7_, provides the dominant contribution. (This is the choice 
that Hartle and Hawking originally made 38 .) Neglecting the prefactor, we therefore 
obtain a no-boundary wavefunction 



* HH (a, 4>) oc exp 



1 



3V(<f>) 



exp 



-1 



W{cj>) 



{l-a 2 V(<P)) 



3/2 



(4.21) 



in the region a 2 V < 1. Using the WKB matching procedure one can show that the 
corresponding solution in the region a 2 V > 1 is 



* (a,0) oc exp 



cos 



1 



WW) (oVw - 1)3/2 - i 



(4.22) 



which is the superposition of the two WKB components of (3.60): 



* = * 

HH - 



± oc exp 



w{4>) 



exp < ±i 



(4.23) 



One may observe that (4.21) and (4.22) agree with the solutions (4.4) and (4.7) 
found earlier by direct examination of the Wheeler-DeWitt equation in the appro- 
priate limits (4.5) and (4.8). Thus although the saddle point corresponding to the 
less than half filled 4-sphere does not appear to be picked out in any special way 
by the path integral, it is favoured by the "approximate" boundary condition 52 ' 119 
that ^ -> 1 as a -> 0. 

One must add the caveat that the approximate boundary condition should be 
amended if one is to consider genuinely complex solutions of the Euclidean field 
equations (4.11)-(4.13). This issue has been considered by Lyons 116 . In general 
one must analytically continue the boundary condition (4.14) demanded by the 
no-boundary proposal. Although the simple picture of matching a real Euclidean 
solution to a real Lorentzian solution at the junction is no longer maintained, one 
nonetheless finds solutions which are initially approximately Euclidean and at late 
times are approximately Lorentzian, with classical inflationary behaviour 116 . 



4-2. The tunneling proposal 

An alternative approach advocated by Vilenkin is that the boundary condition 
for the wavefunction, \I>, should be such as to embody the notion that the universe 
"tunnels into existence from nothing" without making such specific restrictions on 
the "initial" geometry as the Hartle-Hawking proposal does. Conceivably there are 
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many possible ways in which such a notion could be translated mathematically into 
a definition of the wavefunction, and indeed many alternative formulations of the 
tunneling proposal have been put forward^. Some early versions were phrased in 
a similar fashion to the no-boundary proposal: in particular, Vilenkin proposed 
defining the wavefunction by a functional integral over Lorentzian metrics which 
interpolate between a given matter configuration, and 3-geometry, hij, and a 
vanishing 3-geometry, 0, lying to its past 

y[h ij ,$,Z] = J2 DgD^e lS[9 ^^ ] . (4.24) 

M J 

Vilenkin has also given an alternative formulation of the tunneling proposal 
in terms of a boundary condition on superspace 40 ' 41 rather than a restriction on 
manifolds included in the path integral. In order to formulate boundary conditions 
on superspace it is necessary to consider its boundary, which can be thought of as 
consisting of 3-metric and matter configurations for which the 3-curvature is infinite, 
or |$| — > oo etc. As we have already seen from the example of S 4 (Fig. 7), not all 
singular 3-geometries will correspond to singular 4-geometries, as it is possible to 
obtain a singular 3-geometry by a degenerate slicing of the 4-geometry. Therefore 
we should distinguish points on the boundary of superspace which correspond to 
genuine singularities of the 4-geometry from those that correspond to degenerate 
slicings*. We call the former the singular boundary of superspace, and the latter 
the non-singular boundary. 

The tunneling proposal of Vilenkin 41 is that the wavefunction, ^ , should be 
everywhere bounded, and at singular boundaries of superspace $ includes only out- 
going modes, i.e., those that carry a flux out of superspace. Thus ingoing modes 
can only enter at the nonsingular boundary. This definition is somewhat vague as 
there is no obvious rigorous definition of positive and negative frequency modes in 
superspace due to the fact that it possesses no Killing vectors 96 , and thus there 
is no clear notion which modes are "ingoing" and which are "outgoing". Further- 
more, the structure of superspace is not completely understood, and it has not been 
rigorously shown that its boundary can be split into singular and non-singular parts. 

The tunneling proposal has in fact been formulated with the minisuperspace 
WKB approximation in mind, in which case the notion of the boundary of minisu- 
perspace and the notions of ingoing and outgoing modes are more clearly defined. 
Since each oscillatory WKB mode \E'~e ™ has a current (3.44), we can classify 
the modes as ingoing or outgoing according to the direction of VS n on the surface 
in question. Heuristically, the idea underlying the Vilenkin boundary condition is 
that the ensemble of universes described by ^ should not include any universes 



* Linde's proposal [69] embodies a similar philosophy. However, it gives a different wavefunction to 
Vilenkin's "tunneling" wavefunction in simple minisuperspace models [125]. 
* This distinction can be made more precise using Morse functions [42] . 
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contracting down from infinite size 40 , but only those that correspond to "tunneling 
from nothing" . As we shall see, the "outgoing flux" condition 41 accords with this 
notion at least in the case of simple minisuperspace models. 

The outgoing-flux version of the tunneling proposal agrees with the path in- 
tegral formulation (4.24) in the case of the simplest de Sitter minisuperspace ' 66 , 
but the two versions would not appear to be equivalent in general 68 . Whereas the 
no-boundary wavefunction fixes the initial data but leaves the contour of path inte- 
gration ambiguous, the tunneling proposal fixes the contour of integration but leaves 
some ambiguity in the specification of the initial data, even when the outgoing-flux 
condition is imposed 68 ' 126 . 

Consider the minisuperspace model of §3.7. In the diagrams of Fig. 5 all surfaces 
3^ R and the points i° L R and i + will be part of the singular boundary, whereas the 

point i~, which corresponds to a — > (a — > — oo) with finite, is the only point 
of the nonsingular boundary. Of course, the condition that a — > is not enough to 
guarantee regularity of the 4-geometry (4.15), and in general one might not expect 
boundary points to cleanly fall into the category of the "singular boundary" or the 
"nonsingular boundary" . However, as we have already observed in the last section, 
in the present minisuperspace model the additional requirement for regularity that 
g^| = ±1 is guaranteed by the constraint equation (4.11) if a(0) = is imposed. 

The oscillatory WKB region, as shaded in the conformal diagrams of Fig. 5 
and Fig. 9, is always bounded by i° R , !J+ and i + , and in all cases except that of 
Fig. 5(f) v the oscillatory region is bounded by i° L and 3~ also. The wavefunction 
in this region is given by a superposition of terms e lS ™, where S n is a solution to 
the Hamilton- Jacobi equation (3.34), which in terms of the variables (a, </>) may be 
written 

-(I) 2+ (S) 2+u(q ^ - < 4 - 25 > 

where U(a, 0) = 2e 3o W(a, </>) = e 4a (e 2a y(</>) — l), and for convenience in what 
follows we have suppressed the index n. The characteristics of (4.25) satisfy 

da _ d<p _ dS _ d(S, Q ) _ d(S,^) 



2S, a — 2S,^ 2U U, a U,^, 
Thus each S(a,0) describes a congruence of classical paths with 

d</> S,s =fS,0 



_ (4 27) 



Unlike the other cases the potential of Fig. 5(f), which corresponds to a potential in which super- 
symmetry is broken through gaugino condensation in string theory, vanishes as 4> — > — oo. This limit is the 
"weak coupling limit" of string theory. 
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Fig. 9: Some schematic probability flows consistent with the tunneling proposal. 
An indicative oscillatory WKB region is shown (shaded area) for the potential 

v = 4> 2 . 

Since U > in the oscillatory region (assuming V(0) > 0), it follows that these 
integral curves satisfy ^ < 1, i.e., they are "timelike" in the minisuperspace 

coordinates and have an endpoint at i + . Since 7t Q = —^-f^- = S, a it follows that 
the WKB modes correspond to expanding universes (a > 0) with 7T a = S, Q < 0, 
or contracting universes with 7t Q = S, a > 0, if we assume Af > 0. Since all paths 
originate at finite values of a the latter solutions will extend from i + into the interior 
of the minisuperspace, i.e., they are "ingoing modes" and are to be excluded on 
account of the tunneling boundary condition. This boundary condition only allows 
the expanding solutions which are "outgoing" at i + . In the approximation in which 
(4.5) applies it follows that the tunneling proposal demands that only the modes 

3 /2~ 

with \1/ ~exp (e 2a V((j)) — l) are admitted in the oscillatory region. The 

tunneling wavefunction is thus complex in the oscillatory region, in contrast to the 



This choice is a matter of convention and the opposite choice would reverse the roles of the "ingoing" 
and "outgoing" modes. Equivalently, the coordinate t is an arbitrary label from the point of view of general 
relativity, and it is a matter of convention whether we choose t to increase or decrease towards the "future" ; 
the "future" being defined by the expansion of the universe (the cosmological arrow) or the increase of 
entropy (the thermodynamic arrow) . In terms of Fig. 4 it amounts to an arbitrary choice between A + and 
A- as representing the "late-time behaviour" of the inflationary solutions. 
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no-boundary wavefunction (4.22) which is real. 

To extend the tunneling wavefunction into the "tunneling region" (the unshaded 
regions of Fig. 5 and Fig. 9) , we can use the WKB matching procedure. Of course, 
there will also be additional solutions that remain entirely in the tunneling region 
and never cross into the oscillatory region. In the case of the null boundaries, 3~ R , 

which border this region we observe that as a — > — oo (3.59) becomes the wave 
equation in the (a, 4>) coordinates, and thus the solutions are asymptotically null 
which leads to a notion of "ingoing" and "outgoing" modes, so that the tunneling 
condition can be imposed. Some possible probability flows consistent with the 
tunneling proposal are shown in the conformal diagram of Fig. 9. 

5. The Predictions of Quantum Cosmology 

Given the many problems and uncertainties in the quantum cosmology pro- 
gramme that have been discussed above, one could easily form the opinion that it 
is premature to talk about the predictions that quantum cosmology makes. Nev- 
ertheless, it is important to investigate the types of predictions we might expect 
quantum cosmology to make about the universe, as well as realistically evaluating 
the limitations of the these predictions. I will concentrate on three key areas in this 
section. Other topics, most notably the variability of the constants of nature, have 
been discussed recently by Vilenkin 125 ' 127 . 

5.1. The period of inflation 

The construction of a suitable probability measure should allow us to answer 
questions such as whether inflation is a feature of a typical universe. As we saw 
in §3.6, the construction of a general probability measure is problematic. Thus we 
shall limit our discussion to the oscillatory WKB limit in which case the Klein- 
Gordon-type current (3.42) is adequate. 

The issue of the duration of the period of inflation has been a point of some 
debate between proponents of the no-boundary wavefunction and proponents of 
the tunneling wavefunction. Consider the canonical minisuperspace model of §3.7. 
We have seen that in the WKB limit the Hartle-Hawking boundary condition gives 
rise to a wavefunction (4.22) in the oscillatory region, which is strongly peaked 
about the set of classical solutions (3.62) that correspond to an inflating universe: 
a(t) oc e^*, 4>(t) ~ O = const. 

Vilenkin 40-42 ' 70 has also studied the minisuperspace model of §3.7, but with 
a different choice of factor ordering in the Wheeler-DeWitt equation (3.59). He 
concludes that similarly to (4.22) the tunneling boundary condition leads to a WKB 
wavefunction 



\l/ v (a, 4>) oc exp 



( 



3V(4>) 



-1 



) 



exp 



3T(0) 



(aM^)-l) 3/2 + ^ 



(5.1) 
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in the oscillatory region, which is also peaked about the classical trajectories of the 
inflationary universes (3.62). 

Both the no-boundary and tunneling wavefunctions thus predict inflation, at 
least in the context of this simple minisuperspace model. The important question 
is how much inflation do the models predict? This will be largely decided by the 
value, (f> , of the scalar field with which the universe "nucleates" in the semiclassical 
regime. 

In order to study the probability flux arising from (3.42) in the WKB limit 
it is necessary to focus on one particular WKB component. The no-boundary 
wavefunction (4.22) is of course real, and the resulting current (3.42) is identically 
zero. However, \& HH is the superposition (4.23) of two WKB components, which 
will correspond to contracting and expanding universes. As discussed in §3.5 it 
is assumed that a decoherence mechanism exists so that the interference between 
the two components is negligible 89 , and we can therefore assume that the universe 
is peaked about one or other WKB component for the purposes of determining 
the probability measure. The difference between the no-boundary and tunneling 
wavefunctions may therefore not seem great once decoherence to a classical universe 
is assumed. In particular, if we take the "outgoing" WKB component in the no- 
boundary case then the only significant difference between \I/ and \E' T is the <p- 

dependent part of the prefactor, exp (^ 3 y^ ) > which is obtained from boundary 

conditions set in the tunneling region. This difference will have ramifications for 
the probability flux, however. 

One question which we might hope to answer in our minisuperspace model 
would be: given that a Lorentzian universe nucleates, what is the probability that 
it inflates by a sufficient amount (~65 e-folds) to solve the problems of the stan- 
dard cosmology mentioned in §1? The answer to this would involve integrating 
the probability flux (3.45) on the surface separating the tunneling and oscillatory 
regions, which is roughly given by a 2 V((/)) = 1. However, our discussion here is 
limited by the fact that our WKB approximation applies to trajectories with (j> ~ 0. 
For such trajectories the probability current, 3, points chiefly in the direction of the 
a-coordinate in minisuperspace. We can therefore attempt to answer the question 
approximately by evaluating the probability current (3.44) on a surface L in min- 
isuperspace with a = const. For sufficiently large values of a such surfaces will lie 
almost entirely within the oscillatory regime (see Fig. 10). On these surfaces we 
obtain a probability flux 



Although the integral of (5.2) may diverge for particular potentials V((j)) of interest, 
this should not be viewed as a problem since, as was emphasised in §3.6, questions in 
quantum cosmology can only refer to conditional rather than absolute probabilities. 




no-boundary wavefunction, \& 



tunneling wavefunction, 



HH' 



(5.2) 
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In the present case let us assume that inflation occurs for large values of the 
scalar field, as is the case for potentials of the "chaotic" type, V = X(p 2p . There 
will then be a minimum value of the scalar field, <j) s , for which sufficient inflation 
is obtained. For V = X(/) 2p a universe with = O initially will undergo 



e- folds of inflation , so that in the case of the quadratic potential (p = 1) we find 
4> suS > 6.6, for example. The relevant conditional probability for sufficient inflation 
on an a = const surface is then given by 



where the (+) case refers to the no-boundary wavefunction and the (— ) case to the 
tunneling wavefunction, and the values 4> 1 and <p 2 are respectively lower and upper 
cutoffs on the allowed values of <p, which must be determined by physical criteria. 
A minimum cut-off might be expected for a variety of physical reasons, such as 
avoiding classical universes which rapidly recollapse. 

In their original investigations of the question of the duration of the period of 
inflation for the quadratic scalar potential, Hawking and Page 52 took 2 = oo, in 
which case both integrals in (5.4) are dominated by large values of <f>, leading to 
a probability V — 1, and thus a "prediction" of inflation. However, this result 
has been criticised by Vilenkin 41 since for k = +1 we have V = jQirip lanck V so 
that values of m<j) > 4/3 are in excess of the Planck scale, and the semiclassical 
approximation will no longer apply. Vilenkin suggested that an upper cutoff, 2 , 
should be introduced at the Planck scale. If this is the case then provided the 
lower cutoff 4> 1 is sufficiently close to zero, we would find that the integral in the 
denominator of (5.4) becomes very large in the case of the no-boundary wavefunction 
(+ sign), leading to "P <C 1, whereas this would not be the case for the tunneling 
wavefunction (— sign), and thus the latter would predict more inflation. 

Introducing a cutoff at the Planck scale might be deemed a rather arbitrary 
procedure, since without any knowledge of Planck scale physics it is impossible 
to be sure whether the "real" answer is better approximated by the introduction 
of a cutoff or not. However, it has been argued on the basis of investigations at 
the 1-loop level that the wavefunction is damped at large values of <fi by quantum 
corrections 128 ' 129 . This renders the wavefunction normalisable and would justify a 
cutoff, 2 , near the Planck scale. 

^ In the slow-rolling approximation [109] it follows from (3.54) and (3.56) (with JV = 1) that the number 
of e-folds is iV e = 6 d(f> V/V' , where <f> is the initial value and <f> e the final value of the scalar field at the 
end of the inflationary epoch. For V = \(f> 2p we have N e — ^ (0 O 2 — e 2 ) . The value of <j) e can be estimated 
from the limit set by \V /V\ < 6 and \V" /V\ < 9. 




(5.3) 




(5.4) 
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Fig. 10: Conformal diagram for V = 0.04c/) 2 . The oscillatory region, given roughly 
by a 2 V > 1, is lightly shaded. Lines a = const are superimposed. For very large 
values of <j> these lie almost entirely in the oscillatory region. The region of ^-values 
excluded by a Planck scale cutoff is darkly shaded. 

Given that the no-boundary wavefunction apparently yields a wavefunction 
peaked around the lower cutoff, (f) 1 , it is important to determine what a reasonable 
value of this cutoff should be. This issue has been considered by Grishchuk and 
Rozhansky 130 , and also more recently by Lukas 131 , who have conducted numerical 
investigations to analyse the behaviour of the caustic 132 in the (a, <fi) plane which 
separates the Euclidean and Lorentzian solutions. This is illustrated in Fig. 11, 
where classical solutions to the Euclidean field equations (3.54)-(3.56) (with Af = 1) 
are shown. The solutions beginning at a = initially follow lines with <f> ~ 0, but 
the approximation eventually breaks down when the trajectories curl back and rec- 
ollapse with a — > and \(f>\ oc — In a — > oo. They thus represent a flow from i~ to 
J~ in the conformal diagrams. The solutions cross each other on the caustic, which 
for large values of </> corresponds to a 2 V((p) = 1, in accordance with our expectation 
from the WKB approximation. 

To the right of the Euclidean solutions in Fig. 11 we would find Lorentzian 
solutions with <p ~ sufficiently far away from the caustic 132 . These solutions 
are not depicted here. The "nucleation of a universe" would thus correspond to 
a solution of the Wheeler-DeWitt equation which was initially peaked about a 
classical Euclidean trajectory with <f> ~ 0, and which then crossed over the caustic 
to be peaked about a corresponding trajectory with <f> ~ in the Lorentzian region. 
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Fig. 11: Classical Euclidean trajectories for the potential V — </> 4 . The approximate 
caustic a 2 V = 1 is indicated by a dashed line, and the improved caustic by a dotted 
line. (From [131].) 

While crossing the caustic, the wavefunction would be peaked about a complex 
solution which was neither truly Euclidean nor truly Lorentzian. 

It is evident from Fig. 11 that for small values of <p the caustic begins to deviate 
from the curve a 2 V = 1. In fact, one finds 130 ' 131 that for sufficiently small values of 
<fi Lorentzian universes never nucleate at all. Accordingly, a portion of the shaded 
region in the conformal diagrams should be excised about the <p = axis in the 
vicinity of i + . The cutoff, <fi^, at which this occurs has been estimated by Lukas 131 
as being 

h ~ p + [2p [l - (1 + e" 1 )" 1 / 2 ] } V2 , (5.5) 

in the case of the chaotic potentials, V = X(p 2p . Provided that the interpretation 
of (5.4) remains valid*, we should therefore set 4> 1 = 0^ as the lower bound in 
the integrals. In the case of the quadratic potential, V = m 2 (p 2 , we then have 
4> 1 ~ 1.5, and using the values of </> suff and (f> 2 found earlier it is straightforward to 

* As a cautionary note, one should observe that since the approximation of the caustic by the curve 
a 2 V — 1 breaks down for small <f>, the use of the surfaces a — const which was assumed in (5.4) may not be 
appropriate for small <j>, and ideally one should determine the probability flux across the caustic itself. Of 
course, in a more careful analysis one would consider complex solutions to the field equations in line with 
[116], rather than real Euclidean and real Lorentzian solutions with a junction condition. One would hope 
that this would not alter the conclusions of the analysis much. 
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check that for typical values of m <C 1, (5.4) gives V <C 1 for \I/ HH and V ~ 1 for \& v 
in accordance with the earlier discussion. 

The above analysis would appear to indicate that the no-boundary wavefunction 
effectively "predicts" a value of <p ~ <p^ , which in the case of the quadratic potential 
is unfortunately less than 4> suS - However, the calculation is model-dependent, and 
if one could find a potential for which ^ > suff then both \1/ HH and \& v would yield 
"P~l. If we compare (5.3) and (5.5) (for N e ~ 65) we see that in the case of the 
chaotic potentials, V = X(p 2p : this requirement is equivalent to p > 62. However, it 
requires an enormously small value of A to keep such a V below the Planck scale 
and this is not promising. Lukas 131 has also estimated the value of for some 
other potentials but did not find any candidates with 0^ > <^> suff . However, given 
the model-dependence of the calculations it cannot be ruled out that some other 
potential, or the coupled effects of two or more scalar fields, might give ^ > 4> suff 
and thereby make inflation a "prediction" of the no-boundary wavefunction. 

5.2. The origin of density perturbations 

It is expected that the anisotropies in the cosmic microwave background radi- 
ation, which were first definitively observed in 1992, have their origin in quantum 
fluctuations in the very early universe. Such fluctuations can be described using 
the formalism of quantum field theory in curved spacetime. Since the era of quan- 
tum cosmology is in a sense prior to that in which quantum field theory in curved 
spacetime is applicable, one would hope to trace the origin of the primordial per- 
turbations back to quantum cosmology. Indeed, this can be done 133-135 and I will 
very briefly outline the main results. 

With homogeneous minisuperspaces as a starting point, one can add small in- 
homogeneous perturbations to the metric and matter fields: 



Here we have restricted attention to the k = +1 FRW minisuperspace model, the 
subscript zero denotes the unperturbed quantities, and Vtij is the unperturbed stan- 
dard round metric on the 3-sphere. The perturbations can be expanded in terms 
of spherical harmonics on the 3-sphere 133 . If one substitutes the ansatz (5.6)-(5.9) 
into the classical action (2.12) and expands to quadratic order one obtains an action 
which can be split into the original minisuperspace action, S , and an additional 
action, S 2 , quadratic in the perturbations, 



hij (^E, £) 

N{x,t) 

KM 



a 2 (t) (Slij + Eij), 
<f> (t) + 6$(x,t), 

+ SN\(x,t). 



(5.6) 
(5.7) 
(5.8) 
(5.9) 



S = S [g A , Af ] + S 2 [g A , Af , £ij , 5$, 5Af , Mi] , 



(5.10) 
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and a corresponding Hamiltonian 

H=N (h + J d 3 xW 2 j + J d 3 x 5J\fH 1 + J d?x8NiH\ (5.11) 

where H is the unperturbed Hamiltonian (3.54), and the terms H 1 and TC 2 are 
linear and quadratic in the perturbations respectively. There is now a non-trivial 
momentum constraint at each point iGE, H l (x) = 0, while the Hamiltonian con- 
straint splits into a piece linear in the perturbations, r H 1 {x) = 0, plus a homogeneous 
piece 

H + H 2 eH + / d 3 xH 2 = 0. (5.12) 



We may quantise (5.12) in the standard fashion to obtain a modified Wheeler- 
Dewitt equation 

= [-1 V 2 + U(q) + H 2 ] * = 0, (5.13) 

in place of (3.27), where the Laplacian is still defined in terms of the original min- 
isuperspace coordinates according to (3.28), and H 2 is a second order differential 
operator which results from quantisation of the homogeneous part of the perturba- 
tions. In the case of pure scalar field modes, for example, the perturbations may be 
decomposed as 

6*(x,t) = ^J2^lm(t)Q n im(^ (5-14) 
nml 

where the Q n i m satisfy the 3-dimensional Laplace equation on S 3 , 

(3) AQV = -(n 2 - l)Q" lm , (5.15) 

and one finds 1 ' 89 ' 133 

A 2 = \ E (4t^ + [(^ - ^ + ™ V ] f ^n) (5-16) 
nml { tt ° r nlm J 

in the case of the quadratic scalar potential. 

It is possible to find solutions to the Wheeler-DeWitt equation (5.13) in which 
the minisuperspace coordinates, q A , are treated semiclassically in the WKB ap- 
proximation, while the perturbations are treated quantum mechanically. One can 
show 1 ' 89 ' 133 ' 134 that the solutions take the form 

* = C(g,$)e iS o^'*)-ij;, (5.17) 

where S satisfies the unperturbed Hamilton-Jacobi equation (3.34), the prefactor 
C depends only on the unperturbed minisuperspace coordinates, and the functions 
i[) satisfy the functional Schrodinger equation 



= H 2 ^. (5.18) 
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The different modes of the scalar perturbations (5.14) do not interact, for example, 
and in this case 

A> = JJ A>nlm{t, f n l m ), (5.19) 
nml 

where each mode ^ n im separately satisfies (5.17) with H 2 given by (5.16). 

The wavefunction (5.17) is thus peaked about classical trajectories, with correc- 
tions, il>, which satisfy the functional Schrodinger equation along these trajectories. 
This is in fact precisely the starting point for the treatment of matter modes in a 
curved spacetime background using the formalism of quantum field theory in curved 
spacetime. Thus the above result is important in that it demonstrates that quan- 
tum cosmology is consistent with the standard approach to the quantum treatment 
of cosmological perturbations. As an added bonus the imposition of a boundary 
condition on such as the no-boundary or tunneling condition, will result in the 
choice of particular solutions of the functional Schrodinger equation, and conse- 
quently a particular vacuum state for the matter modes. Both the no-boundary 
condition and the tunneling condition pick out 41 ' 136-139 a de-Sitter invariant state 
known as the "Euclidean" or Bunch-Davies vacuum 140 ' 141 , which is the state that 
is often assumed in cosmological calculations of density perturbations. In fact, this 
state is picked out by many boundary conditions 137 ' 138 , and thus could be regarded 
as a natural quantum state for matter in quantum cosmology. 

One may solve the functional Schrodinger equation and study the growth 133 of 
the modes. As discussed above, the results obtained are the same as those which 
are found if one begins with quantum matter fields in de-Sitter invariant vacuum 
states in background inflationary spacetimes 142 ' 143 . In particular, one obtains a 
scale-free spectrum of density perturbations which act as seeds for the formation 
of galaxies and other structures in the universe. Such a spectrum accords well 
with the spectrum deduced from the COBE measurements of cosmic microwave 
anisotropies 144 . 

5.3. The arrow of time 

The question of the origin of the arrow of time in the face of the CPT-invariance 
of the laws of physics is one of the deepest unresolved conceptual issues in theoretical 
physics, and it has been the subject of more than one major conference 4 ' . The 
question of the nature of time assumes prime importance in quantum gravity and 
quantum cosmology^. One might argue that if the basic laws of physics are time 
symmetric, then the arrow of time could have its origin in a boundary condition for 
the wavefunction of the universe. 

There are many different arrows of time observed in the universe. Some of 
these such as the psychological arrow of time (we remember only the past), and 

The conceptual issues surrounding the nature of time in quantum cosmology would entail a series of 
lectures in themselves. For further discussion see, e.g., [24], [146]-[148] and references therein. 
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the electromagnetic arrow of time (the choice of retarded as opposed to advanced 
solutions of Maxwell's equations), could well be argued to be consequences of the 
thermodynamic arrow of time, which arises from the second law of thermodynamics. 
However, the expansion of the universe provides an alternate cosmological arrow of 
time which is not obviously directly related to the increase of entropy. Whether 
such a relationship does exist is a question which might potentially be resolved by 
quantum cosmology. 

In 1985 Hawking 149 proposed that the thermodynamic arrow and cosmological 
arrows of time were correlated: that is to say if the universe were spatially closed 
then entropy would decrease in the contracting phase®. In arriving at this proposal 
he was influenced by the fact that the no-boundary wavefunction is CPT-invariant, 
and also by some early studies of simple minisuperspace models which possessed 
quasi-periodic solutions 151 which "bounced" instead of recontracting to a singular- 
ity as a — > 0. However, Hawking soon changed his mind on the issue, which he 
terms his "greatest mistake" 152 . This change of mind was brought about by a num- 
ber of factors. Firstly, Page pointed out that CPT-invariance of the no-boundary 
wavefunction does not imply CPT-invariance for an individual WKB component of 
the wavefunction, which would correspond to the history of a classical universe 153 . 
Furthermore, minisuperspace models which were subsequently studied, such as that 
of the Kantowski- Sachs universe, were found not to admit bounce solutions but al- 
ways possessed singularities to the future 154 ' 155 . Finally, as was mentioned in §4.1, 
the bounce solutions do not feature even in simple minisuperspace models once 
one considers the contribution of complex metrics 116 . In general, the approximate 
minisuperspace boundary condition that ^ — > 1 as a — > must be altered to al- 
low for approximately Euclidean metrics which also contribute a rapidly oscillating 
component to the wavefunction 116 ' 156 . 

Hawking, Laflamme and Lyons 156 have recently argued that a thermodynamic 
arrow of time results from the imposition of the Hartle-Hawking boundary condi- 
tion. More precisely, they have considered the evolution of primordial fluctuations 
as outlined in §5.2, but accounting for the "approximately Euclidean" geometries 116 
which appear to be required if one considers complex metrics. They find 156 that 
the gravitational wave perturbations have an amplitude that remains in the linear 
regime and is roughly time-symmetric about the time of maximum expansion. Such 
perturbations cannot be said to give rise to an arrow of time. Density perturba- 
tions behave differently, however. They start out small but grow large and become 
non-linear as the universe expands, and moreover this growth continues during the 
contracting phase of the spatially closed universe. This growth of inhomogeneity 
therefore provides an arrow of time which could be considered to be an essentially 
thermodynamic arrow. Since it does not match the cosmological arrow in the con- 
tracting phase, the only reason for the coincidence of the two arrows in the present 
epoch would appear to be an anthropic one. In particular, the conditions that would 

® This idea was first suggested by Gold [150], and has resurfaced a number of times in different contexts. 
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prevail in a contracting phase would appear to preclude the existence of life 152 ' 156 . 
Thus the fact that we are around to make observations means we must find ourselves 
in a cosmological epoch in which the two arrows coincide. 

It should be added that the debate about whether the cosmological and ther- 
modynamic arrows of time coincide has not yet been closed, however, and it is 
still maintained by Kiefer and Zeh 157 that the boundary condition on the Wheeler- 
DeWitt equation must be such that the thermodynamic arrow would reverse in a 
recontracting universe. 



6. Conclusion 

I hope to have shown you that although research in quantum cosmology is still 
rather speculative and open-ended, its framework nonetheless has the potential not 
only to provide answers to questions surrounding the origin and early evolution of 
the universe, but also to help us unravel the mysteries of quantum gravity. Quantum 
cosmology is a field in which a great deal remains to be done, and the results of 
§5 should be regarded as a hopeful indication of the types of predictions we might 
hope to make. It is too early to draw definitive conclusions about the relative merits 
of the various boundary condition proposals. The results of §5.1 ostensibly favour 
the "tunneling proposal" . However, this conclusion is only based on a few simple 
models, and therefore some caution must be exercised. 

As a final note, it is worth mentioning that recent results show that super- 
symmetry provides a means of restricting possible boundary conditions for the 
Wheeler-DeWitt equation, or the corresponding Dirac square-root equation. The 
supergravity constraint equations for various homogeneous minisuper space models 
appear to be so restrictive that they only pick out the most symmetric quantum 
states 102 ' 103 ' 158 . In particular, simple analytic solutions for \1/ in the supersymmet- 
ric Bianchi-IX minisuperspace have been found 102 in the empty and filled fermion 
sectors, which have a natural interpretation as 158 wormhole states 121 ' 159 , or as 160 
Hartle-Hawking no-boundary states. Some doubts were initially cast on the rele- 
vance of these solutions, as the states appear to have no counterpart in 4-dimensional 
supergravity 161 ' 162 . However, more recent work 46 on supersymmetric minisuper- 
spaces corresponding to Bianchi "class A" 76 models indicates that infinitely many 
physical states with finite (even) fermion number can be found, and these are direct 
analogues of physical states in full supergravity. This result shows that such minisu- 
perspace models are likely to be physically very important for quantum cosmology. 
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